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ABSTRACT 

We present the results of a spectroscopic survey of the kinematic structure of star-forming galaxies 
at redshift z ~ 2 — 3 using Keck/OSIRIS integral field spectroscopy. Our sample is comprised of 
12 galaxies between redshifts z ~ 2.0 and 2.5 and one galaxy at z ~ 3.3 which are well detected in 
either Ha or [O in] emission. These galaxies are generally representative of the mean stellar mass 
of star forming galaxies at similar redshifts, although they tend to have star formation rate surface 
densities slightly higher than the mean. These observations were obtained in conjunction with the 
Keck laser guide star adaptive optics system, with a typical angular resolution after spatial smoothing 
~ 0.15" (approximately 1 kpc at the redshift of the target sample). At most five of these 13 galaxies 
have spatially resolved velocity gradients consistent with rotation while the remaining galaxies have 
relatively featureless or irregular velocity fields. All of our galaxies show local velocity dispersions ~ 60 
- 100 km s" 1 , suggesting that (particularly for those galaxies with featureless velocity fields) rotation 
about a preferred axis may not be the dominant mechanism of physical support. While some galaxies 
show evidence for major mergers such evidence is unrelated to the kinematics of individual components 
(one of our strongest merger candidates also exhibits unambiguous rotational structure), refuting a 
simple bimodal disk/merger classification scheme. We discuss these data in light of complementary 
surveys and extant UV-IR spectroscopy and photometry, concluding that the dynamical importance 
of cold gas may be the primary factor governing the observed kinematics of z ~ 2 galaxies. We 
conclude by speculating on the importance of mechanisms for accreting low angular-momentum gas 
and the early formation of quasi-spheroidal systems in the young universe. 

Subject headings: galaxies: high-redshift - galaxies: kinematics and dynamics - galaxies: starburst 
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1. INTRODUCTION 

Galaxies are the discrete luminous building blocks of 
the visible universe, tracing the development of gravita- 
tional structures across cosmic ages. The earliest galax- 
ies likely formed near the density peaks of the primor- 
dial power spectrum (e.g., Bardeen et al. 1986) which 
were able to decouple from the cosmic expansion at early 
times and form a burst of stars from the primordial gas. 
As time progressed however, these chaotic proto-galaxies 
gradually evolved, merging with their neighbors in newly 
collapsing dark matter halos, accreting greater quantities 
of gaseous fuel from a filamentary intergalactic medium, 
and polluting their environments with the metallic detri- 
tus of their early stellar generations. It is at intermediate 
redshifts z ~ 2 — 3 that these morphologically irregular, 
juvenile galaxies are thought to have formed the major- 
ity of the stellar mass which we observe in modern-day 
galaxies (Dickinson et al. 2003, Reddy et al. 2008). 
Through a combination of galaxy-galaxy mergers, rapid 
star formation, and secular evolution, these galaxies ex- 
perienced a strong morphological transformation into the 
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coherent structures of the Hubble sequence which have 
predominated since redshift z ~ 1 (Giavalisco et al. 1996; 
Papovich et al. 2005). 

Despite our growing knowledge of the broad global 
characteristics of high redshift galaxies (e.g., Erb et al. 
2006c; Reddy et al. 2006a; Papovich et al. 2006; Cowie 
et al. 2008; and references therein), our understanding of 
their internal structure and dynamical evolution has been 
limited by their small angular size, typically < 1 arcsec- 
ond. Such objects are generally not well resolved by the 
ground-based imaging and spectroscopy which form the 
backbone of our observational data. This limitation pre- 
cludes us from addressing such questions as: (i) what 
are the triggering and regulation mechanisms of their 
starbursts? (ii) does star formation occur in the (cir- 
cum) nuclear regions of dynamically relaxed systems or as 
a result of tidal shocks induced by major mergers? (iii) 
do individual star-forming regions follow a global abun- 
dance pattern or are there strong variations in chemical 
enrichment within a galaxy? Each of these distinctions 
has implications for the evolution and development of 
structure and stellar populations within a given galaxy. 

One of the most promising means to investigate these 
questions is studying the ionized gas surrounding bright 
star-forming regions produced by the flood of energetic 
photons from young stars. Excited atoms in this gas lose 
energy predominantly by emission of strong rest-frame 
optical nebular emission lines (redshifted into the NIR 
for galaxies at z ~ 2 — 3) such as Ha, H/3, [On], [Oni], 
and [Nil]. The relative strengths of these lines encode 
information about the chemical composition of the gas 
and the shape of the ionizing spectrum, permitting de- 
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duction of the metallicity and star formation rates of the 
emission regions. These emission features also have nar- 
row natural line widths and are good kinematic tracers 
of the ionized gas. 

Such efforts using slit spectroscopy (e.g., Pettini et al. 
2001 at z ~ 3; Erb et al. 2004, 2006b, 2006c at z ~ 2; 
Weiner et al. 2006 at z ~ 1) have suggested that the 
kinematics of these galaxies may frequently be inconsis- 
tent with simple rotationally supported gas disk models. 
However, these studies have been limited by the small an- 
gular size of typical galaxies relative to the atmospheric 
seeing. Often there may be only 1-2 spatially indepen- 
dent samples across the face of a given galaxy, and ad- 
ditional uncertainty may be introduced by misalignment 
of the slit with the (a priori) unknown kinematic axis. It 
is therefore unclear whether the high observed velocity 
dispersion is genuine or caused by the smearing of unre- 
solved kinematic structure, such as disk-like rotation or 
merging clumps. 

The recent advent of adaptive optics (AO) on 10 m 
class telescopes has offered the opportunity to overcome 
the limitations previously imposed by atmospheric tur- 
bulence by rapidly correcting the distorted wavefront us- 
ing deformable mirrors. Paired with integral-field spec- 
troscopy (IFS), it is possible to obtain diagnostic spectra 
of spatial regions resolved on scales of order 100 mil- 
liarcseconds (mas), corresponding to roughly 1 kpc at 
rcdshift z ~ 2 — 3. The data provide an empirical answer 
to whether the velocity fields of galaxies in the early uni- 
verse are predominantly represented by virialized disk- 
like systems, major mergers, or some other dynamical 
structure, and whether the resulting star formation is 
uniform in its properties across a given galaxy or exhibits 
regional variation on kiloparsec scales. 

Early results from such IFS programs with (Genzel et 
al. 2006; Law et al. 2007a; Wright et al. 2007, 2009) 
and without (Forster Schreiber et al. 2006; Bouche et 
al. 2007; Genzel et al. 2008; Nesvadba et al. 2008; van 
Starkenburg et al. 2008) 7 adaptive optics, in addition 
to a small number of galaxies observed with the aid of 
gravitational lensing (Nesvadba et al. 2006; Swinbank 
et al. 2007; Stark et al. 2008; Jones et al. in prep.), 
have helped confirm that the high velocity dispersion is 
an intrinsic property of high-redshift galaxies. The total 
number of individual galaxies observed however is ex- 
tremely small, and considerable disagreement remains as 
to the underlying nature of the dynamical structure of 
these galaxies, with Genzel et al. (2008; and references 
therein) suggesting that rotationally supported gaseous 
disks form the majority of the bright galaxy population 
while Law et al. (2007a) favor a significantly lower rota- 
tion fraction with the majority of galaxies kinematically 
dominated by their high velocity dispersions. 

In this paper we greatly expand the original sample 
of Law et al. (2007a) by presenting spatially resolved 
laser guide star adaptive optics (LGSAO) spectroscopy 
of an additional 10 galaxies at redshifts 2.0 < z < 2.5. 
With the larger sample size it is possible to draw broader 
conclusions about the kinematical characteristics of the 
galaxy population. The paper is structured as follows. In 

7 One galaxy presented by Forster Schreiber et al. (2006) was 
observed with the aid of adaptive optics, as were two of the five 
main galaxies discussed by Genzel et al. (2008). 



$2] we describe our sample selection, observational tech- 
nique, and data reduction algorithms. In fj3]we discuss 
the physical properties of our galaxies, particularly their 
morphological f£ )3.ip and kinematic f £)3.3[) structure and 
stellar populations ( §3.2| . The characteristics of each of 
our 13 galaxies are discussed individually in detail in |21 
fJ5] discusses the characteristics of the global population, 
comparing to recent observations by other groups in £15 . 21 
and exploring possible implications for gas accretion and 
formation mechanisms in H5.3l and l5.4l We conclude with 
a brief summary in Sj6j while a derivation of the statistical 
properties of an ensemble of inclined systems is included 
in an appendix ( §A.lj) . We assume a standard ACDM 
cosmology based on 3-year WMAP data (Spergel et al. 
2007) in which H Q = 73.2 km s" 1 Mpc" 1 , Q m = 0.238, 
and n A = 0.762. 

2. OBSERVATIONS 
2.1. Target Selection 

Recent years have witnessed an abundance of meth- 
ods for locating galaxies in the redshift range z ~ 2 — 3. 
These methods include optical (U n GlZ) color selection 
(e.g., Steidel et al. 2003, 2004), near-IR BzK (Daddi 
et al. 2004) and J - K (Franx et al. 2003) color se- 
lection, selection by sub-mm flux density (Chapman et 
al. 2005), and Lya emission surveys (e.g., Martin et al. 
2008; and references therein). Of these samples, those 
arising from the optical color selection are perhaps the 
most well studied and account for most of the global star 
formation activity at z ~ 2 — 3 (Reddy et al. 2005). 
As described in detail by Adelberger et al. (2004), the 
U n G1Z selection technique represents a generalized ver- 
sion of the Lyman-break technique employed by Steidel 
et al. (2003) to identify rapidly star-forming galaxies at 
z ~ 3 on the basis of their strong 912 A Lyman break, 
rcdshifted into the U n bandpass. In this paper, we fo- 
cus primarily on this U n GlZ sample, and particularly on 
those galaxies in the redshift range z = 1.8 — 2.6 (i.e., the 
"BX" galaxy sample of Steidel et al. 2004). While initial 
identification of these galaxies was based upon photomet- 
ric preselection, extensive rest-UV spectroscopic follow- 
up work has derived precise redshifts for all galaxies in 
our target sample. 

The galaxies selected in this manner are typically 
bright and actively forming stars, with mean extinction- 
corrected star formation rates (SFR) ~ 30 Af yr _1 (see 
discussion by Erb et al. 2006b), and SFR surface densi- 
ties similar to those observed in local starburst galaxies 
(e.g., Kennicutt et al. 1998b). The resulting winds from 
supernovae and massive stars drive energetic (~ a few 
hundred km s _1 ; Steidel et al. in prep.) large-scale 
outflows into the IGM surrounding these galaxies, cre- 
ating the ubiquitous blueshifted interstellar absorption 
features observed in rest-frame UV spectra (e.g., Pettini 
et al. 2002; Shapley et al. 2003), and the corresponding 
redshifts of the resonantly scattered Lya emission line 
(Verhamme et al. 2008). With the aid of deep near-IR 
(e.g., Erb et al. 2006c) and mid-IR (e.g., Shapley et al. 
2005b; Papovich et al. 2006; Reddy et al. 2006b; and 
references therein) photometry, stellar population mod- 
eling suggests that galaxies at z ~ 2 — 3 span a broad 
range of stellar masses and evolutionary states. 

Individual galaxies selected for study were drawn from 
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this pool of available targets subject to a variety of cri- 
teria. Some targets were deliberately selected for their 
young stellar population ages and correspondingly small 
stellar masses, some for their old ages and large stel- 
lar masses, and some for other reasons including com- 
plex or multi-component rest-UV morphologies (Law et 
al. 2007b), strong detections in Ha narrowband surveys 
(Q1700-BX710 and Q1700-BX763), unusual spectral fea- 
tures, or previous acquisition of long-slit kinematic data. 
Given the relatively shallow OSIRIS X-band sensitiv- 
ity limit (SFRi im ~ 1 M Q yr" 1 kpc~ 2 in 2 hours of in- 
tegration; Law et al. 2007a) however, the most com- 
mon criterion applied was preferential selection of those 
galaxies for which previous long-slit spectroscopy (Erb 
et al. 2006b) indicated nebular line fluxes > 5 x 10~ 17 
erg s" 1 cm -2 . We discuss the nature of our final target 
sample further in i J5.ll 

Additional physical constraints require that there be 
a suitably bright (R < 17) star within ~ 60" of each 
galaxy to use as a tip-tilt (TT) reference for the LGSAO 
system 8 , and that the wavelength of redshifted line emis- 
sion falls between the strong spectroscopically unresolved 
night-sky OH emission features that dominate the near- 
IR background, and avoids wavelengths of extremely 
strong telluric absorption bands. The first of these con- 
straints will not systematically bias the resulting sample 
given its reliance only on the distribution of nearby stars, 
while the second imposes blackout "windows" in redshift 
space, the most significant of which spans z ~ 2.6 — 2.9. 
The final target sample is listed in Table Q] and in- 
cludes 24 galaxies in 9 distinct fields distributed widely 
across the sky. Of the 24 galaxies observed, 13 galax- 
ies (distributed amongst 7 fields) have significant detec- 
tions while an additional two (Q1623-BX455 and Q2343- 
BX587) are too poorly detected to merit inclusion in our 
analysis. 

2.2. Observational Technique 

Observations were performed using the OSIRIS 
(Larkin et al. 2006) integral-field spectrograph in com- 
bination with the Keck II LGSAO system during 6 ob- 
serving runs between June 2006 and September 2008. 
The majority of these sessions followed the basic observ- 
ing scheme described in §2.1 of Law et al. (2007a), to 
which we refer the reader for a detailed description of 
our approach. In brief, we observed the target galax- 
ies in narrowband filters (bandwidth ~ 0.1 (im) corre- 
sponding to the wavelength of either Ha or [O in] A5007 
(hereafter [O in] ) emission lines using the 50 mas lenslet 
sampling scale with 900 second individual exposures. 
Given the redshift range of most of our target galax- 
ies (z ~ 2 — 2.5), Ha is redshifted into the if-band 
and we generally choose to observe this emission line 
since the quality of the LGSAO correction is highest at 
longer wavelengths. At wavelengths longer than ~ 2.26 
fun, however, the thermal continuum background seen by 

The theoretical sky coverage of 

the Keck LGSAO system fsee, .e.g., 

http:/ /www2. keck. hawaii.edu/optics/lgsao/performance. html ) 
is ~ 60%, ranging from almost 100% near the Galactic plane to 
~ 20% or less near the Galactic poles and in particularly sparse 
fields such as the Hubble Deep Field. Generally, we found that ~ 
70 - 80 % of our desired targets were close enough to a bright star 
to be observed with the LGSAO system. 



OSIRIS becomes comparable to that of the atmospheric 
OH lines, significantly degrading the quality of the ob- 
servational data. We therefore observe [O in] emission 
in either the H or K bands as necessary for some galax- 
ies (e.g., Q1623-BX543 and DSF2237a-C2). The obser- 
vational configuration for each of our target galaxies is 
detailed in Table [1] 

During a typical observing sequence we first acquire 
a brief (~ 60 second) observation of our TT reference 
star in order to center our pointing and provide both a 
point-spread function (PSF) and flux calibration refer- 
ence. When possible, we obtained another observation 
of the TT reference star at the end of our integration 
on a given galaxy in order to bracket the science obser- 
vations and measure any changes in the PSF. We offset 
from the TT star to the target galaxy using precise off- 
sets measured from deep (~ 10 hour) ground-based op- 
tical and/or HST Advanced Camera for Surveys (ACS) 
imaging data (where available). Each galaxy was ob- 
served in pairs of 900 second exposures dithered within 
the field of view, typically by ~ 1.4" along the long axis. 
This dithering maximizes our on-source integration time 
but effectively halves our useful field of view to ~ 1" x 
1.5" (still significantly larger than the typical galaxy, see 
discussion by Law et al. 2007b). Each exposure pair 
was repeated with small ~ 50 mas dithers for a total of 
between 1800 and 12600 seconds of integration. Total 
integration times for each object (Table [T]) were dictated 
by observational availability and the need to obtain suf- 
ficient integration to ensure high-quality detections. It 
was typically possible to confirm detection of a given 
galaxy in the difference of two 900 second exposures, al- 
lowing us to determine in real-time the targets which 
were likely to produce the best results given sufficient 
integration. With some notable exceptions (e.g., Q2343- 
BX389) deeper integrations were not performed if a given 
galaxy was undetected within 1-2 hours. 

2.3. Data Reduction and Flux Calibration 

Data reduction was performed using a combination of 
OSIRIS pipeline and custom IDL routines described in 
§2.2 of Law et al. (2007a) which produce a compos- 
ite, three-dimensional data cube (consisting of two image 
dimensions and one spectral dimension) for each target 
galaxy. Given the multi-year baseline over which our ob- 
servations were obtained, the reduction algorithms have 
evolved over time to reflect the changing performance 
characteristics of OSIRIS. The most significant modifi- 
cation was introduced for data obtained in June 2007 
(i.e., for our observation of galaxies HDF-BX1564 and 
Q1623-BX502) which were affected by an imperfect fo- 
cus within OSIRIS resulting in a PSF slightly elongated 
along an axis roughly 45 degrees to the lenslet grid. This 
was corrected by implementing a CLEAN-type algorithm 
(e.g., Hogbom et al. 1974) to deconvolve the elongated 
PSF and replace it with a more typical circular PSF. 

As described by Law et al. (2007a) the composite data 
cube for each galaxy is sub-sampled by a factor of two 
in each spatial dimension and convolved with a Gaus- 
sian kernel in order to increase the signal-to-noise ratio 
(SNR) of the spectrum in each spatial pixel (spaxel). The 
width of this Gaussian is chosen as required to produce 
the highest quality smoothed data cube without signif- 
icantly inflating the PSF delivered by the LGSAO sys- 
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tern. We typically chose a FWHM of 80 mas, although 
for some extremely low surface brightness galaxies (e.g., 
HDF-BX1564) we increased this smoothing to ~ 200 mas 
in order to increase the SNR at the cost of spatial resolu- 
tion. The effective width of the PSF (as measured from 
the TT star) before and after this smoothing is given for 
all galaxies in Table ffl For one galaxy (HDF-BX1564) 
we also smooth the spectrum of each lenslet in the wave- 
length domain by a kernel of FWHM 5 A to better dis- 
tinguish the faint emission feature from the surrounding 
noise. We note that while these algorithms are largely 
automated, the details for each observation are carefully 
optimized to produce the highest-possible quality final 
products. 

Raw composite spectra of each galaxy (shown in Fig. 
[T]) were obtained by summing into a single spectrum the 
data from all spaxels in a box encompassing the nebular 
emission morphology, thereby sampling both bright and 
faint emission regions. This raw spectrum is calibrated 
using a telluric transmission spectrum of the night sky 
to normalize the throughput as a function of wavelength. 
This calibration spectrum is determined directly from ob- 
servations of bright telluric standard stars for data ob- 
tained in June/September 2008, and from the OSIRIS + 
Keck LGSAO system -I- atmospheric models of Law et 
al. (2006) for all previous data (see discussion in Law 
et al. 2007a). The absolute flux calibration of the cor- 
rected spectra is determined by matching extracted spec- 
tra of the TT stars to IR photometry given in the 2MASS 
Point-Source-Catalog. We note that this method differs 
from that adopted by Law et al. (2007a), and slight dif- 
ferences from their calibration are to be expected. 

The largest single uncertainty in the flux calibration 
of these data arises from the nature of the LGSAO PSF, 
which can vary both with isoplanatic angle and (rapidly) 
with time, resulting in significant fluctuations in the per- 
centage of total light residing in the AO-corrected core 
of the PSF. These fluctuations are poorly understood, 
and correlated primarily with mid-level atmospheric tur- 
bulence not well corrected by the AO system. In ad- 
dition low surface-brightness or extremely broad spec- 
tral features may be missed entirely (see discussion in 
§5.ip . leading to an underestimate of the total flux from 
galaxies with appreciable flux in such regions. Given the 
uncertainties in our bootstrapped flux calibration, we es- 
timate that the systematic flux uncertainty for any given 
source is ~ 30% (see also discussion in Law et al. 2007a). 

In Table [2] we list the [Oin], Ha, and [Nil] emission 
line fluxes of each galaxy (with uncertainties based on 
the noise present in the underlying spectra after sub- 
traction of Gaussian models for the identifiable emission 
components 9 ) along with estimates of the oxygen abun- 
dance based on the Pettini & Pagel (2004) calibration 
of the N2 — Fryq u i /-Fka relation. Most of these fluxes 

are about half the aperture-corrected values determined 
from previous long-slit NIRSPEC spectroscopy (Erb et 
al. 2006b). One possible explanation for this discrep- 
ancy may be that the factor of 2 aperture correction used 
to calibrate the NIRSPEC data could overestimate slit 
losses (particularly for the most compact galaxies), re- 

9 Residual features at the emission line wavelength (after sub- 
traction of the Gaussian model profile) are generally indistinguish- 
able from the noise in the rest of the spectrum. 



suiting in erroneously high flux estimates. Alternatively, 
this may also be due to low surface brightness emission 
which OSIRIS is not sufficiently sensitive to detect. 

Total nebular emission line luminosity is calculated 
from these fluxes assuming a standard cosmological 
model and correcting for extinction derived from stellar 
population models f H3.2[) using a Calzetti et al. (2000) 
attenuation law modified as described by Erb et al. 
(2006c) 10 . We convert this extinction-corrected luminos- 
ity to an Ha star-formation rate using the Kennicutt et 
al. (1994) calibration: 

SFR(M yr-^ i 26x L i ^ rgs _ 1 x 0.56. (1) 

where the factor of 0.56 converts to the Chabrier (2003) 
IMF. For those two galaxies observed in [O in] instead of 

Ha we either assume that . Lh ° = 1 (Q1623-BX543, 

L [Om] 

for which the actual flux ratio is unknown) or adopt 

, Lh " = 0.52 (DSF2237a-C2, see discussion in Law et 
L [Om] 

al. 2007a). We also list the central wavelength of the pri- 
mary observed emission feature in each spectrum. Sys- 
temic redshifts are determined by dividing these wave- 
lengths by the rest-frame vacuum wavelengths of Ha and 
[Om] (i.e., 6564.614 A and 5008.239 A, respectively) 
and correcting for the heliocentric motion of the Earth 
at the date, time, and direction of observation (calcu- 
lated using the IRAF rvcorrect algorithm). We do not 
include fainter [N n] or [O in] A4960 in our systemic red- 
shift calculations as these features are generally too weak 
to improve the redshift fit. 

Kinematic data were extracted from the oversampled 
composite data cubes using custom IDL routines to fit 
Gaussian profiles to the spectrum at each spatial loca- 
tion, thereby obtaining maps of the emission line flux, 
wavelength centroid, and spectral FWHM across each 
galaxy. This information is converted to maps of the ve- 
locity relative to the systemic redshift and the velocity 
dispersion (with the instrumental resolution [R ~ 3600] 
subtracted off in quadrature) as shown in Figure O As 
discussed by Law et al. (2007a), the quality of these 
maps is improved considerably by the spatial smooth- 
ing described above, and by excluding any fits which are 
either unphysical or for which the data had inadequate 
SNR. 

3. RESULTS 

3.1. Morphologies 

As demonstrated by numerous authors (e.g. Abraham 
et al. 1996; Conselice et al. 2005; Elmegreen et al. 2005; 
Lotz et al. 2006; Law et al. 2007b) the morphological 
structure of galaxies at redshift z > 2 is highly irregular 
and frequently composed of multiple spatially separated 
components when observed in the rest-frame UV. It is 
worthwhile exploring the similarity between such rest- 
UV morphologies and those of the ionized gas. Given 
the extremely faint nature of these galaxies, we produce 
maps of the ionized gas morphology (Fig. [21 left-hand 
panels) by summing spectral channels within lcr of the 
peak emission wavelength at each oversampled spatial 

10 E(B - V) ranges from 0.005 to 0.370. 
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Fig. 1. — OSIRIS spectra integrated over the spatial extent of each galaxy. Vertical dashed lines indicate the fiducial locations of nebular 
emission lines based on the systemic redshift. Hashed green regions indicate the locations of strong atmospheric OH emission features 
which can give rise to strong residuals in the spectra. 

sure of the number of components of the light distribu- 
tion, see Law et al. 2007b), the second-order moment 
M 2 o of those pixels constituting the top 20% of the total 
flux (Lotz et al. 2004), the total luminous area / and ef- 
fective radius r corrected for the PSF as described in Law 
et al. (2007a), and the distance d,2 C between individual 
components (for galaxies which have spatially separated 
pieces). These measurements are given in Table [3J 

We note, however, that it is difficult to compare these 
values directly to rest-UV morphological data presented 
by Lotz et al. (2004, 2006) and Law et al. (2007b) due 
to systematic differences in the observational data. First 



pixel meeting the SNR criteria described in M3.31 Such 
a 'tunable' narrow-band filter for each position in the 
galaxy minimizes the contribution of noise to the flux 
map . 

Following the method outlined by Abraham et al. 
(2007) we define a segmentation map for each galaxy 
using a quasi-Petrosian isophotal cut (77 = 0.2) which 
is independent of morphology. Using this segmentation 
map we calculate a variety of morphological statistics in- 
cluding the Gini coefficient G (a measure of the curve of 
growth of the light distribution, see discussion by Abra- 
ham et al. 2003), the multiplicity parameter '5 (a mea- 
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(and most important) the quasi-Petrosian isophotal cut 
simply selects all pixels above our S /N ratio threshold for 
nearly all of our galaxies since our method of construct- 
ing the lowest-noise flux maps has artificially eliminated 
all of the fainter galaxy pixels for which it was not possi- 
ble to fit a reliable emission-line spectrum. While this is 
clearly not ideal from a morphological standpoint, for our 
lowest surface-brightness sources (e.g. HDF-BX1564) 
this spectral line-fitting method and associated SNR cut 
is frequently the only mechanism by which we can reli- 
ably distinguish any features from the background noise 
at all. Our isophotal cut generally does not reach sur- 
face brightness levels as low as those of many narrow- 
band studies; this results in lower values of G than are 
normally derived in the rest-frame UV (e.g., Lotz et al. 
2004, 2006; Law et al. 2007b) as there is less contrast 
among pixels in the segmentation map. In contrast, the 
multiplicity parameter \& is relatively robust to such sur- 
face brightness variations and still reliably indicates the 
presence of multiple components in the light profile. Ad- 
ditionally, many of the galaxies observed are only a few 
times the size of the observational PSF (which can vary 
considerably from galaxy to galaxy) which also has an 
impact on the numerical classification of the morphol- 
ogy. We therefore recommend that these morphological 
statistics be interpreted in a relative, rather than abso- 
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lute, sense. 

Four of our target galaxies (HDF-BX1564, Q1700- 
BX490, Q1700-BX710, and Q1700-BX763) lie in fields 
for which deep optical imaging data has been obtained 
with HST-ACS. These data in the HDF and Q1700 fields 
have been described by Law et al. (2007b) and Peter et 
al. (2007), respectively. The comparative morphologies 
of HST-ACS (i.e., tracing the rest-UV continuum) and 
OSIRIS Ha emission are shown in Figure [3] Given the 
factor of ~ 5 difference between the limiting SFR sur- 
face density probed by the HST-ACS and OSIRIS data 
(~ 0.2M Q yr^ 1 kpc~ 2 and ~ lAf© yr~ x kpc~ 2 respec- 
tively; see discussion in Law et al. 2007a) there is gener- 
ally a good correspondence between the respective mor- 
phologies. Due to the extremely narrow field of view of 
OSIRIS, there are no absolute reference points to which 
we can calibrate the coordinate solution as our offsets 
from the TT acquisition star are only reliable to within 
~ 100 mas (i.e., comparable to our PSF) due to various 
uncertainties resulting from (for example) measuring the 
coordinates of the TT star, global uncertainties in the 
coordinate solution of the ground based imaging data, 
and the unknown proper motion of the TT stars. It is 
not possible therefore to reliably align the HST-ACS and 
OSIRIS images to high precision using absolute coordi- 
nates alone. Instead, we align the HST-ACS contours 
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to target as needed to contain the emission line regions. The FWHM of the PSF (after smoothing described in i|2,3H is indicated in the 
left-hand panel for each galaxy. Contours represent linear intervals in line flux density. All images are presented in a standard orientation 
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shown in Figure [3] with the OSIRIS maps "by-eye", slid- 
ing the contours freely until they appear to best over- 
lap the Ha flux data. While this prohibits us from in- 
vestigating small discrepancies between the locations of 
peak emission in isolated single sources, we are still able 
to compare spatial differences for sources with multiple 
emission features and the relative shapes of the various 
emission profiles. We discuss the specific features of in- 
dividual galaxies in H4. II ~ |4~T31 below. 

3.2. Stellar Population Modelling and Gas Masses 

Using our extensive ground-based U n GlZJK s and 
Spitzer IRAC and MIPS photometry in the selected sur- 
vey fields it is possible to construct stellar population 
models for the spectral energy distributions of the target 
galaxies. The modeling procedure is described in detail 
by Shapley et al. (2005a). In brief, we use Bruzual and 
Chariot (2003) models with a constant star formation 
rate (unless otherwise specified), solar metallicity, and a 
Chabrier (2003) initial mass function (IMF). These mod- 
els are overplotted on our photometric data in Figure [4] 
and the individual values for SFR, stellar mass (M*), 
population age, and extinction are tabulated in Table 
[5l As discussed by previous authors (e.g., Shapley et al. 
2001, 2005a; Papovich et al. 2001) , degeneracies in the 
stellar population model tend to cancel out with regard 
to the total stellar mass, meaning that this is typically 
the best constrained of these parameters. We note that 
previous modeling of many of these galaxies has been pre- 
sented by Erb et al. (2006c) and Shapley et al. (2001); 
the previous models are updated here for consistency and 
to add Spitzer IRAC data where available. 

The galaxies for which the poorest fit to the photom- 
etry is obtained (HDF-BX1564, Q1623-BX502, Q1623- 
BX543, Q2343-BX660) tend to be those with multiple 
components to the light distribution. In general, how- 
ever, there is reasonable agreement between the values 
of SFR derived from SED fitting and those implied by 
the Ha luminosity (Table [5]) , although the SED-based 
estimates tend to be systematically greater than those 
based on Ha. 

We also include in Table [5] estimates of the mass in 
cold gas of each of our galaxies (and the corresponding 
gas fraction /i = ^/ gas /(M gas + M*)) determined by Erb 
et al. (2006c) based on long-slit NIRSPEC spectroscopy. 
While it would be possible to use the global Schmidt law 
(Kennicutt 1998) to calculate new estimates of M gas from 
our observed Ha fluxes, the results of this calculation are 
sensitive to both Ha flux and the apparent size of the 
galaxy, which systematically differ between OSIRIS and 
NIRSPEC due to the different surface brightness thresh- 
holds reached by the two spectrographs. We use the Erb 
et al. (2006c) estimates for consistency in our later ( §5.2| 
comparison to other galaxy samples. 

3.3. Kinematics 

As illustrated by Figure [2] (middle panels), galaxies 
show a mix of velocity fields ranging from largely fea- 
tureless (Q1623-BX453) to smoothly varying (Q1700- 
BX490) to extremely complex (Q2343-BX513). In par- 
ticular, we find that five of our galaxies (Q1623-BX502, 
Q1623-BX543, Q1700-BX490, DFS2237a-C2, Q2343- 
BX418) have velocity fields consistent with rotation- 
like velocity gradients, six (Q0449-BX93, Q1217-BX95, 



HDF-BX1564, Q1623-BX453, Q1700-BX710, Q1700- 
BX763) consistent with negligible velocity structure, 
and two (Q2343-BX513, Q2343-BX660) with velocity 
structure unlike simple rotation. Additionally, three of 
our 13 galaxies (HDF-BX1564, Q1623-BX543, Q1700- 
BX490) consist of multiple spatially (and kinemati- 
cally) distinct regions: two of these galaxies display 
significant kinematic shear in their individual compo- 
nents (Q1623-BX543, Q1700-BX490) while one does not 
(HDF-BX1564). 

We discuss the kinematics of individual galaxies in de- 
tail in §4]below. The single feature common to all of these 
galaxies however is the significant velocity component 
with no preferred kinematic axis (Fig. [21 right-hand pan- 
els) which dominates over all coherent kinematic struc- 
tures. This velocity dispersion may be quantified in two 
ways. An estimate of the overall velocity dispersion cr ne t 
of a galaxy may be obtained by simply fitting a single 
Gaussian profile to the spatially integrated spectrum of 
the entire galaxy (i.e., the spectra shown in Fig. [1]). This 
measurement does not distinguish between spatially re- 
solved velocity gradients and small-scale motions within 
a given resolution element. While <7 ne t can provide a 
reasonable estimate of the dynamical mass using the for- 
mula 

Af dyn = ( 2 ) 

(where C — 5 for a uniform sphere; see Erb et al. 2006c), 
it is not an optimal means of determining the intrin- 
sic dispersion of the ionized gas. We therefore also cal- 
culate cr m0 an which is the flux- weighted mean 11 of the 
velocity dispersions measured in each individual spaxel 
(i.e., the 'local' velocity dispersion shown in the right- 
hand panels of Figure [2]) ■ This quantity effectively sup- 
presses spatially resolved velocity gradients and provides 
a more accurate measure of the typical line of sight ve- 
locity dispersion at a given location within the galaxy, 
albeit with some potential bias due to the flux weight- 
ing, beam smearing, and finite sampling of the observa- 
tional data. Values of tx me an and <r ne t are tabulated in 
Table [H In the three cases where the galaxy consists of 
spatially distinct regions we calculate values for cr m can in 
each of the regions and also measure the kinematic off- 
set V2c between them. As expected, values of cr mcan are 
generally consistent with or slightly less than a DCt except 
for cases where large-scale velocity structure inflates a net 
(e.g. HDF-BX1564). 

We place these descriptive arguments on a firmer nu- 
merical footing by constructing one-dimensional velocity 
curves of each of our target galaxies as shown in Figure 
03 These are obtained by extracting spectra from our cal- 
ibrated data cubes along simulated slits whose width is 
matched to the spatial resolution of the data, and orien- 
tations (indicated by red lines in Fig. [2]) chosen to max- 
imize the apparent velocity gradient across the galaxy. 
We define the shear velocity w s hcar as half the maximum 
difference between any two positions along these pseudo- 
slits, noting that this will disregard whether the maxi- 
mum difference occurs between locations at the extreme 
ends of a galaxy or indeed whether the resulting velocity 

11 We adopt a flux -weighted mean rather than a simple mean in 
order to minimize contamination from the lowest S/N ratio spectra. 
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curves resembles actual rotation in any way. Any appar- 
ent trends can easily be confirmed by inspection of the 
two-dimensional velocity map (Fig. [2]) . 

Ideally, it would be possible to perform detailed kine- 
metric analyses of each of our galaxies (e.g. Shapiro et 
al. 2008; Genzel et al. 2008; Wright et al. 2009) in 
order to obtain a set of best-fitting models for the light 
distribution and ionized gas kinematics, along with the 
dynamical ratio Wdisk/cint of the model (where i>disk is 
the maximum circular speed attained within the galaxy 
and Cint is the intrinsic local velocity dispersion). How- 



ever, the observational data generally do not justify such 
detailed analyses because the majority of the galaxies do 
not show concrete evidence (such as spatially resolved 
velocity shear) to suggest that they might be described 
by traditional inclined-disk models. Of the 5 galaxies 
with smoothly varying velocity gradients consistent with 
rotation, only two have these gradients aligned with the 
morphological major axis (as might often be expected for 
a foreshortened disk). Of these two, one (Q1623-BX543) 
has a close (~ 7 kpc) companion which does not partake 
in the same pattern of velocity shear. A model for the 
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single remaining galaxy, DSF2237a-C2, which is isolated 
and has a significant kinematic gradient aligned with its 
morphological major axis has been presented by Law et 
al. (2007a). 

Absent such models, the best kinematic statistic avail- 
able is the observed ratio i>shear/cm ea n, va l ues °f which 
range from 0.2 — 0.7 with a typical uncertainty ~ 0.2 (Ta- 
bled]). While we have attempted to maximize the value 
of ^shear/cmoanby appropriate choice of the position angle 
for our pseudo-slit, and by using the smaller of our two 
estimates of a, the meaning of any particular value must 
be considered with regard to several important factors. 

1. The observed velocity curves will be foreshortened 
by a factor of sin i. On average this will introduce 
a correction factor of n/A (see lA.l[) . 

2. The OSIRIS data are relatively shallow and probe 
radii < 2 kpc. If the velocity curve rises at larger 
radii, our observations will underestimate the true 
circular velocity of the system. 

3. Unresolved kinematic structure smeared by the ob- 



servational beam will decrease the observed veloc- 
ity gradient and inflate the central velocity disper- 
sion to which the flux-weighted <7 mcan is most sen- 
sitive. In the case of a rotating disk model this ef- 
fect will be particularly pronounced at small radii 
where the velocity gradient is steepest. 

All of these effects will systematically produce values 
of Wshoar/o-moan lower than v disk /a int , although the mag- 
nitude of these effects is difficult to quantify as they de- 
pend strongly upon the (unknown) intrinsic kinematic 
structure of the galaxies. Rotating disk models may 
expect correction factors ^2 — 4 (the disk model for 
DSF2237a-C2 developed by Law et al. 2007a, for ex- 
ample, had a correction factor ~ 2), depending on the 
size of the observational PSF and the precise distribu- 
tion of flux with respect to the underlying rotation field. 
Even allowing for an extreme correction factor ~ 5, the 
observed range Vshear/omean = 0.2 — 0.7 is significantly 
lower than the typical values (t>disk/omt ^ 10 — 20, i.e., 
"shcar/cmcan^ 2 — 4) appropriate to disk galaxies today 
(e.g. Dib et al. 2006), in agreement with the general con- 
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elusion that these galaxies must be considerably thicker 
than traditional disk galaxies in the local universe (see 
also discussion by Forster Schreiber et al. 2006, Genzel 
et al. 2006, Law et al. 2007a, Genzel et al. 2008). 

Such corrections to the observational data may not ob- 
viously be desirable however for galaxies which are not 
well-described by inclined disk models. In at least six 
of our 13 galaxies both the velocity and velocity disper- 
sion are consistent with a constant value across multiple 
resolution elements. In the absence of coherent veloc- 
ity structure, the high velocity dispersions observed for 
z ~ 2 galaxies ((<r mcan ) — 78 ± 17 km s _1 ) cannot be 
explained by the convolution of such structure with the 
observational PSF 12 (as, for instance, in the z ~ 1.6 sam- 
ple of Wright et al. 2009), and instead indicate that both 
the high velocity dispersion and low bulk motions about 
a preferred kinematic axis may be intrinsic properties of 
the ionized gas. 

4. NOTES ON INDIVIDUAL GALAXIES 

4.1. Q0449-BX93 

This intermediate-mass galaxy (JVf* = 0.9 x 1O 1O M0) 
is dominated by a single emission component at the sys- 
temic redshift with a faint secondary component offset by 
~ 150 km s _1 located ~ 1 kpc to the northeast. This faint 
secondary feature may represent a kinematically distinct 
star-forming region or a small galaxy in the process of 
merging with the the brighter system (see discussion by 
Law et al. 2007a). Weak [Nil] emission (both A6549 
and A6585 lines are detected with similar significance) is 
roughly concentric with the Ha emission and implies a 
metallicity 12 + log(0/H) = 8.37 ± 0.09 (approximately 
half-solar). There is no evidence for significant bulk mo- 
tions with Vshcar/c mean — 

0.2 ±0.2. 

12 Numerical modelling suggests that beam-smearing can ac- 
count for at most a few percent of the observed (cr mcan ). 



4.2. Q1217-BX95 

This galaxy is qualitatively similar to Q0449-BX93, 
and is well-detected but compact (only slightly larger 
than the observational PSF) with no evidence of ex- 
tended structure despite its moderate stellar mass (3 x 
10 10 Mq). There is no obvious velocity structure within 
this source. 

4.3. HDF-BX1564 

HDF-BX1564 is among the least well detected of our 
galaxies, and required heavy spatial and spectral smooth- 
ing of the composite data cube in order to successfully 
discern emission features. This extensive processing has 
the unfortunate consequence of increasing our effective 
PSF to ~ 0.3", comparable in size to the two primary 
emission regions (northwest [NW] and southeast [SE]) 
comprising the source. These two features are sepa- 
rated by «2c ~ 171 km s -1 and d,2 C — 7 kpc. While a 
third component is apparent in Figure [2 based on anal- 
ysis of typical noise patterns within the composite data 
cube we believe that this feature is probably not real. 
Supporting this explanation is the rest-UV morphology, 
which we show for comparison in Figure [3] Both the NW 
and SE features have counterparts in the ACS imaging 
data, while the questionable northeastern feature does 
not. The HST-ACS imaging data also suggests faint dif- 
fuse emission in the field of view, particularly in the form 
of spurs reaching between the bright knots. Given that 
the bright knots themselves are only barely detected (at a 
S/N ratio of ~ 5) in the heavily smoothed OSIRIS data, 
it is not surprising that we do not detect these spurs, 
even if they have Ha counterparts. 

The local velocity dispersions of the NW and SE fea- 
tures are 60 ± 13 and 86 ± 8 km s _1 respectively, consis- 
tent with the typical values observed for the rest of the 
galaxy sample, while their effective PSF-corrected radii 
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( §3.ip are 1.3 and 0.1 kpc (i.e., the SE feature is con- 
sistent with a point source given the size of the PSF). 
The respective dynamical masses of the two components 
are 0.5 and 0.1 x 10 10 A/ Q . Given the poorly resolved na- 
ture of the SE component however, the radius of this ob- 
ject is particularly ill-defined and given the similar values 
of cr mcan for the two components they may be approxi- 
mately equal mass. 

The multiple components and low surface-brightness 
features in the UV morphology may indicate that HDF- 
BX1564 is an (~ equal mass) major merger (although 
see discussion by Law et al. 2007b) In this case, addi- 
tional insight may be gained from Keck/NIRSPEC long- 
slit spectroscopy obtained in June 2004 (see Erb et al. 
2006b) with a 0.76" wide slit aligned along the axis be- 
tween the two major components to within 5°. As illus- 
trated in FigurelHl both primary components are well de- 
tected in the spectrum with similar kinematic and spatial 



separation as in the OSIRIS data. In this considerably 
deeper observation there is evidence for nebular emission 
extending between the two components, the SE feature 
has a spatial FWHM of 1.1" and reaches all of the way to 
the NW component, whose 0.6" FWHM is more closely 
matched to the PSF of these seeing-limited data. No- 
tably, there is no gradual shift in wavelength along the 
axis connecting the two, but instead a sharp transition 
upon reaching the NW feature, suggesting that the emit- 
ting material between the two components does not par- 
take in a smoothly varying velocity curve and bolstering 
the case for a major merger interpretation of this galaxy. 

Intriguingly, there is faint detection of [Nil] emission 
(12+ \og{0/H) = 8.67 ± 0.04) in the OSIRIS data cen- 
tered on neither of the two primary Ha emitting regions, 
but on a region of space between them corresponding to 
a faint spur in the rest-UV morphology (Fig. [3]). This 
may suggest the presence of an obscured AGN at the 
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dynamical center of this system (although neither the 
rest-UV spectrum nor the SED shows evidence for an ac- 
tive nucleus), while the UV and nebular line emission are 
dominated by off-center regions of relatively unobscured 
star formation. Alternatively, this could be evidence for 
diffuse shocked gas, or low ionization-parameter gas as 
discussed by Martin et al. (1997). Given the low SNR 
of our data however, it is not possible to conclusively 
explain the observational complexities of this system. 

4.4. Q1623-BX453 

As discussed in Law et al. (2007a) this galaxy is one of 
the best examples of the "typical" kinematics observed 
in our sample, with spatially well-resolved nebular emis- 
sion exhibiting a negligible velocity differential across the 
majority of the galaxy and kinematically dominated by 
a high local velocity dispersion cx mcan = 78 ± 23 km s _1 . 
The ^shear/cmean value calculated for this galaxy is dom- 
inated by a region of high relative velocity at the north- 
eastern end of the galaxy. This region is detected with 
low confidence however (SNR ~ 4) and represents an 
abrupt kinematic discontinuity rather than an extended 
gradient. 

Q1623-BX453 is inferred to have a large amount of 
cold gas, a sizeable stellar mass (~ 3 x 1O 1O M0), high 
metallicity (12+ \og[0/H) = 8.67 ± 0.02), and a high 
SFR surface density with a correspondingly large outflow 
velocity ~ 900 km s _1 as traced by rest-UV spectroscopy 
of interstellar absorption lines. Given the large quantity 
of cold gas in this galaxy, it is particularly intriguing that 
there is no evidence for the rotational structure in which 
such gas is traditionally supposed to reside. 

4.5. Q1623-BX502 

Q1623-BX502 is one of the lowest stellar mass objects 
observed with OSIRIS (M* = 3 x lO 9 Af ) and is the only 
galaxy to be successfully detected under poor observing 
conditions with high humidity and V-band seeing > 1" . 
These sub-optimal observing conditions are reflected in 
the relatively broad PSF (220 mas) compared to the 



other galaxies in Figure [2j In addition, Q1623-BX502 
was one of two galaxies which required correction for a 
spatially distorted PSF as described in fc|2.3l While there 
is evidence for resolved velocity structure with peak am- 
plitude ~ 50 km s _1 , this shear is not aligned with any 
obvious morphological major axis. A similar result was 
derived from the SINFONI observations of this galaxy by 
Forster Schreiber et al. (2006). 

4.6. Q1623-BX543 

This galaxy consists of two spatially distinct compo- 
nents separated by «2 C = 125 km s _1 and c?2c = 6.7 kpc 
(0.8") in projection (with one almost due north of the 
other). Both components are present in our deep ground- 
based 7?.-band imaging, which shows a slight elonga- 
tion along the N - S axis, indicating that both compo- 
nents correspond to rest-UV continuum emission regions. 
Given the typical ~ 1" PSF of this rest-UV imaging data 
however it is not possible to determine whether the con- 
tinuum image consists of a single elongated feature or 
two separate features similar to those shown in Figured! 

While not apparent in Figure [2] due to the limitations 
imposed by the dynamic range of the velocity scale, there 
is velocity structure within the northernmost (N) com- 
ponent. In Figure [7] we plot an expanded view of each of 
the components which better shows this structure. While 
the general trend is along a roughly NW-SE axis (see also 
Fig. [5]) there are appreciable deviations from a smooth 
trend along this axis. We note particularly that the rela- 
tive velocity of the S component is strongly inconsistent 
with a smooth continuation of this velocity field (i.e. it 
would need to be blueshifted relative to systemic rather 
than redshifted). Given these discrepancies it is clearly 
impossible to fit both components of this system with a 
single simple kinematic model. 

Perhaps the most likely scenario is that the primary 
(N) component is akin to the typical single galaxies ob- 
served in the rest of our sample, which in this case is ex- 
periencing a merger with the smaller S component. The 
N/S components have <7 me an = 139 and 60 km s^ 1 , and 
effective radii r = 1.1 and 0.7 kpc respectively. Assum- 
ing that each component is relaxed and kinematically 
dominated by <7 mcan we obtain (Eqn. [2]) estimates of the 
dynamical masses Md yn ~ 2.5 and 0.3 x 10 10 M©. If this 
system is a merger, these rough calculations suggest that 
the mass ratio of the merger is ~ 8/1. 

While it is somewhat curious that the photometry of 
this source is (relatively) well-fit by a simple stellar pop- 
ulation model, such a mass ratio suggests that the larger 
of the two components may constitute the bulk of both 
the current and past star formation in the system. The 
total stellar mass of the system is quite small (4x 10 9 M Q ) 
with an extremely young stellar population as might be 
expected if rapid star formation 13 is occurring as a result 
of the rapid injection of large quantities of cold gas. 

4.7. Q1700-BX490 

As in the case of Q1623-BX543, this galaxy has a com- 
plex, multiple-component structure. As demonstrated by 

13 The high SED-derived SFR given in Table[5]for Q1623-BX543 
and Q1700-BX490 are likely significantly overestimated, and the 
10 Myr ages highly uncertain due to difficulties fitting appropriate 
E(B — V) for these very young galaxies (see discussion in Reddy 
et al. 2006b). 
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Fig. 7. — The velocity map of Q1623-BX543 (left panel, see also Fig. |2j is expanded about each of the two components labelled North 
('N') and South ('S'). The color scale represents line-of-sight velocity relative to the systemic redshift in all three panels, the limits of this 
velocity scale vary as indicated to best illustrate the velocity features of each component. 



Figure [31 the overall location and size of these two com- 
ponents is broadly consistent with HST-ACS rest-UV 
imaging. At the fainter surface brightness levels probed 
by HST-ACS both components are enclosed in a com- 
mon envelope of flux, and it is probable that the discon- 
tinuous Ha morphology is simply a consequence of our 
limiting surface brightness sensitivity. Curiously, the W 
component has a single well-defined peak to the UV flux 
distribution located between the double peaks of nebular 
Ha emission (see contours in Fig. [3]). 

In the larger western (W) component there is a 
smoothly varying velocity gradient along the NE - SW 
axis (i.e. not aligned with the E-W major axis of the 
HST-ACS morphology) from which emission line flux ap- 
pears to be concentrated in two closely separated knots. 
Rather than peaking in the center of the object as might 
be expected for a rotating gaseous disk (although see dis- 
cussion by Wright et al. 2009), the velocity dispersion in- 
stead peaks in one of these two knots. While the smaller 
eastern (E) component is close to the main body in pro- 
jection (3.4 kpc) and would likely appear to connect with 
the W component in deeper spectroscopic data, it does 
not partake in the velocity gradient defined by the W 
body. Rather, it has its own (albeit weak) velocity gra- 
dient aligned nearly perpendicular (~ 70°) to that of the 
main body. 

Paired with the additional faint structure evident in 
the HST-ACS image (note particularly another low 
surface-brightness component SW of the main body) one 
natural explanation for this galaxy may be that it repre- 
sents a merger (mass ratio ~ 5/1) with multiple sites 
of star formation. The broadband photometry is ex- 
tremely well fit by a young (10 Myr) stellar population 
model with a small stellar mass (4 x 10 9 Mq) and a high 
current star formation rate. This galaxy is also likely 
quite metal-poor as [Nil] is undetected, placing a limit 
12 + log(0/H) < 8.26 (< 2/5 solar). 

4.8. Q1700-BX710 

As for Q1623-BX453, this galaxy is relatively well- 
resolved with no strong evidence for spatially resolved 
velocity structure. While [Nil] is detected in emission, 
it is weak (12 + \og(0/H) = 8.40 ±0.08) and consistent 
with a point source roughly concentric with the Ha flux 
peak. There is a slightly redshifted extension of flux to 
the northwest of the galaxy core which has a lower ve- 
locity dispersion and is coincident with the low surface 



brightness emission tail seen in HST-ACS imaging data 
(Fig. |5J). We note that another galaxy (Q1700-BX711) is 
located - 5" (41 kpc) to the northwest of Q1700-BX710 
(i.e., approximately along the direction of the 'tail' indi- 
cated by the HST-ACS data) and at the same redshift to 
within 10 km s _1 . However, a long-slit NIRSPEC spec- 
trum obtained along the line connecting the two galaxies 
shows no evidence for star formation occuring in the re- 
gion between the two galaxies. If these galaxies are gen- 
uinely close in space as implied by their rcdshifts they 
may merge over the next few hundred Myr. 

4.9. Q1700-BX763 

Q1700-BX763 has an apparent western spur to the Ha 
flux consistent with that observed in HST-ACS imaging 
data. While there is mild evidence for a velocity gradient 
along this spur the magnitude of this gradient is small 
and its significance low (Fig. |5J). As for Q1700-BX710, 
weak [N n] emission is detected (implying 12 + \og(0/H) 
= 8.41±0.25) that is consistent with a point source and is 
roughly (within ~ 0.1") concentric with the Ha emission 
peak. 

4.10. DSF2237a-C2 

DSF2237a-C2 is located at a significantly higher red- 
shift than any other galaxies in the target sample (z ~ 
3.3), and is drawn from the "LBG" galaxy catalog of Stei- 
del et al. (2003). Given its high redshift, it is interest- 
ing that DSF2237a-C2 has some of the strongest velocity 
shear observed in our sample ( Ushear/omean = 0.6 ± 0.2) 
and is the only isolated case where this shear is consis- 
tent with rotation and unambiguously aligned with the 
morphological major axis. Nonetheless, the local velocity 
dispersion of DSF2237a-C2 is comparable to that of the 
rest of the galaxy sample with tJ mcan = 89 ± 20 km s _1 . 
We note that this galaxy was discussed previously by 
Law et al. (2007a), although various estimate of the 
kinematics and flux parameters have been updated here 
for consistency with our new calibration routines. 14 

4.11. Q2343-BX418 

14 In particular, the value of £rQ m ] in Table 3 of Law et al. 

(2007a) was incorrectly given as 2.2 X 10 42 erg s — 1 , while it should 
instead have been given as 22 X 10 42 erg s — 1 . The revised number 
is consistent with the estimate of 16.8 X 10 42 erg s — 1 derived using 
our new flux calibration routine. 
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Fig. 8. — Flux map of Q2343-BX513 overlaid with contours 
showing the distribution of flux in the blue and red sides of He* 
emission (green and red contours respectively). The wavelength 
difference between the two sets of contours is 10 A, or 150 km s — 1 at 
the redshift of the galaxy. The black cross indicates the location 
of peak [Nil] emission. 

As for Q1217-BX95, Q2343-BX418 is well detected in 
the composite OSIRIS data cube but is spatially compact 
with a slight elongation along the N-S axis. While there 
is some evidence for coherent velocity structure aligned 
with the morphological minor axis, the magnitude of this 
shear is significantly smaller than the local velocity dis- 
persion. We note additionally that this is our lowest 
mass galaxy with M* = 1 x 1O 9 M , a similarly small 
dynamical mass estimate M^ yn = 3 x 10 9 M Q and very 
little dust (E(B - V) = 0.03). With so little extinction, 
Q2343-BX418 is also a strong Lya emitting source (see 
Pettini et al. in prep.). 

4.12. Q2343-BX513 

Q2343-BX513 is elongated along the N-S axis and ex- 
hibits a complex velocity structure composed of multi- 
ple spatially superimposed components. At first glance 
the one-dimensional velocity curve in Figure [5] (with 
^shear/cmean = 0.7 ± 0.2) gives an indication of smooth 
resolved velocity shear in this source, but the two- 
dimensional velocity map (see Figure [2]) is inconsistent 
with a simple rotational model. Following the velocity 
pattern from North to South, we find that the north- 
ernmost region is blueshifted by —80 km s _1 (relative to 
the systemic redshift); moving South we reach a maxi- 
mum redshift of v = +50 km s _1 near the flux peak, but 
then the velocity falls off again to lower values, down to 
v ~ km s _1 at the southern edge of the galaxy. Indeed, 
if we consider the spatial distribution of the red and blue 
sides of the Ha emission individually we note that while 
the redder component is well-centered the blue edge of 
Ha emission is dominated by two distinct emission re- 
gions located directly north and south of the center (Fig. 
I!]). The total Ha flux of Q2343-BX513 is dominated by 
the centrally concentrated redshifted component. 

This double-component morphology helps to explain 
the previous long-slit observations of Erb et al. (2006c) 
who found two slightly different redshifts for the system 
with different a when the galaxy was observed at two 
different position angles. Notably, the galaxy also has a 
highly unusual double-featured UV spectrum with Lya 
emission at redshifts z — 2.106 and 2.114 (i.e. offset 



by 770 km s _1 ) and possible doubles of the strongest 
interstellar absorption lines, such as Si II A1260, Oi + 
SnA1302, and CnA1334. This spectrum is additionally 
confused by the presence of absorption lines arising in 
an intervening system Q2343-MD80 (z = 2.014) located 
a few arcseconds away. [N n] A6583 is detected with a 
strength relative to Ha which implies 12 + log(0/H) = 
8.63 ± 0.03, close to the solar oxygen abundance, and is 
roughly concentric with the redder Ha component with 
a similar north-south elongation. There is a significant 
velocity gradient in the [Nil] emission with an end-to- 
end velocity differential ^110 km s _1 . 

We note that this galaxy is the oldest in the OSIRIS 
sample (among both detected and undetected targets) 
with an estimated stellar population age of ~ 3 Gyr and 
a large stellar mass M* = 7.8 x 10 10 M©. One possibility, 
therefore may be that the redshifted component (with 
associated [N n] emission and the larger velocity disper- 
sion) is the massive "mature" galaxy which contains the 
majority of the evolved stellar population and has an ap- 
preciable velocity gradient (as expected for such a mas- 
sive galaxy; see discussion in Sj5|. Behind this component 
may be a smaller (infalling) system which is partially ob- 
scured by the central core of the foreground galaxy, caus- 
ing the apparent two-component morphology. While it 
is beyond the scope of the present work to verify this 
hypothesis (and indeed, the two components may simply 
be unrelated juxtapositions in redshift space), we may 
certainly conclude that this system is particularly com- 
plex. 

4.13. Q2343-BX660 

The Ha morphology of Q2343-BX660 is composed of 
two spatially distinct peaks (which are more apparent in 
the 3-d data cube than in the collapsed Ha map) which 
blend together at low surface brightness to form a con- 
tiguous system. While there is a strong velocity differen- 
tial across the combined system, the line-of-sight veloc- 
ity remains relatively constant across each of these two 
clumps and changes rapidly in the region between the 
two (~ 80 km s _1 over less than 0.1"). While turnovers 
in the velocity profile are naturally found in rotating sys- 
tems beyond a certain radius, it is suspicious that the en- 
tire velocity gradient of this galaxy is contained within 
a single angular resolution element (corresponding to a 
radius ~ 400 pc). 

Instead, this figure is reminiscent of Figure 8 in Law 
et al. (2006), whose panel (d) shows the nearly identi- 
cal velocity map expected for a 2-clump velocity model 
in which the apparent scale of the 'turnover' is deter- 
mined by the observational PSF blurring together emis- 
sion from the two regions. Indeed, Nesvadba et al. (2008; 
see also Wright et al. 2009) recently presented observa- 
tions of a galaxy at z ~ 3 (Q0347-383 C5) which showed 
a very similar velocity structure, transitioning rapidly 
from a plateau of ~ 60 km s _1 to ~ —80 km s _1 within 
a negligible physical distance, leading the authors to fa- 
vor a merger interpretation for the galaxy. Similarly, we 
conclude that the kinematics of Q2343-BX660 are more 
consistent with merging compact galaxies than with un- 
derlying rotation in a single relaxed system. Given the 
relative constancy of a across the galaxy we surmise that 
in such an event the mass ratio should be < 3/1, with a 
total stellar mass of ~ 1O 1O M0. 
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Fig. 9. — Mean He* surface brightness as a function of stellar 
mass for the Erb et al. (2006c) long-slit survey of galaxy kinemat- 
ics. Data are plotted for all galaxies in the Erb et al. (2006c) survey 
in the redshift range z = 1.9 — 2.5 which have well-defined values of 
both the Ha flux and galaxy size. Open circles represent galaxies 
which have been successfully observed with OSIRIS, crosses galax- 
ies which were either undetected or poorly detected with OSIRIS, 
and filled squares galaxies not observed with OSIRIS. 

5. DISCUSSION 

5.1. Selection Effects 

As outlined in §2.11 targets were selected from the 
rest-UV galaxy sample for a wide variety of reasons 
(high/low stellar mass, extended/compact morphology, 
presence/absence of kinematic shear in long-slit spectra, 
etc.) subject to the general criterion that we believed the 
galaxy would be detected on the basis of previous long- 
slit observations (Erb et al. 2006b) or narrowband Ha 
imaging. While our original choice of targets attempted 
to include a wide range of galaxy properties, it is still dif- 
ficult to know how representative our results are of the 
general population of star-forming galaxies at z = 2 — 3, 
given the large fraction of targets which were undetected 
or poorly detected (11 out of 24) despite extensive ob- 
servational efforts (see Table [lj . This bias is difficult to 
quantify since the probability of detecting a given galaxy 
is a strong function of observing conditions; 11 of the 16 
unsuccessful observations (some of which represent mul- 
tiple attempts to observe the same galaxy) listed in Table 
Q]were performed under sub-optimal conditions typically 
characterized by poor seeing > 1" (resulting in poorer 
AO correction), high humidity (exacerbating telluric ab- 
sorption), and/or attenuation by cirrus (> 1 mag). In 
contrast, only one of our successfully detected galaxies 
(Q1623-BX502) was observed under similarly poor con- 
ditions. 

The likelihood of detection was probably also governed 
by the Ha surface brightness of the target galaxies. In 
Figure [9] we plot the surface brightness £ = -Fka/^Ha °f 
galaxies targeted by the Erb et al. (2006c; based on val- 
ues given in their Table 4) long-slit spectroscopic survey, 
indicating galaxies which were detected/ undetected in 
the OSIRIS observations. While the results naturally will 
be quite inhomogenous given the widely varying weather 
conditions and exposure times of the OSIRIS observa- 



tions we note that the mean surface brightness of galax- 
ies in the NIRSPEC survey which were successfully de- 
tected with OSIRIS is 4.3 x 10~ 18 erg s" 1 cm" 2 kpc~ 2 , 
compared to 2.3 x 10 -18 erg s _1 cm -2 kpc~ 2 for galaxies 
either undetected or poorly detected with OSIRIS, and a 
mean/median of 2.6/1.7 xl0~ 18 erg s _1 cm" 2 kpc~ 2 re- 
spectively for the entire NIRSPEC sample. Clearly, we 
are observing a population of galaxies with somewhat 
higher than average star formation rate surface density. 

Such a surface-brightness limit will also bias our ob- 
servations towards high surface-brightness regions within 
a given galaxy. Indeed, evidence suggests that at least 
some of the observations may genuinely be missing some 
flux from fainter surface-brightness regions. First, the in- 
tegrated Ha fluxes are about half the aperture-corrected 
values determined from long-slit spectroscopy by Erb et 
al. (2006b), although in some cases this may also indi- 
cate that the aperture correction adopted for the long- 
slit data was overestimated. Similarly, galaxy sizes es- 
timated from the long-slit program tend to be some- 
what larger (although considerably less certain given the 
seeing-limitcd PSF) than those from OSIRIS. Addition- 
ally, the dynamical mass estimates are considerably less 
than the total of the stellar + gas masses, as might be 
expected if the radius of the system is underestimated 
in Eqn. [2] (or the stellar and/or gas masses are dis- 
tributed over a larger volume than the ionized gas traced 
by OSIRIS). 

In some cases (e.g. DSF2237a-C2) the kinematics sug- 
gest that the rotation curve may continue to rise at larger 
radii, meaning that deeper spectroscopy could result in 
larger values for u s hcar- This is unlikely to be the case 
in general however, since at least at least six of the 13 
galaxies are consistent with a velocity gradient of = 
km s _1 kpc -1 . Indeed, deeper imaging of one galaxy 
(Q1623-BX502) with VLT/SINFONI (Forster Schreiber 
et al. 2006) showed no evidence for significant flux or 
velocity shear beyond that seen with OSIRIS. 

Figure [9] also suggests that we are more likely to 
detect younger galaxies than older galaxies with a 
larger established stellar mass. The 13 galaxies de- 
tected with OSIRIS have a mean assembled stellar 
mass (log(M <t /M )} = 10.1, compared to the mean 
( log (M*/M©)) = 10.6 (i.e., three times more massive) 
of the eleven galaxies which were undetected (or only 
poorly detected) in the OSIRIS data. Thus, while our 
sample may be representative of typical-mass galaxies 
in the star- forming population (the full sample of 818 
optically-selected galaxies at redshifts 1.8 < z < 2.6 with 
stellar mass information in our various survey fields has 
a mean stellar mass (log(M»/M )) = 10.07; Fig. \W± see 
also Reddy et al. 2006a), the observations presented here 
are not particularly sensitive to the high-mass end of the 
star-forming galaxy distribution. We note for complete- 
ness however that all of these masses are high relative 
to the total galaxy population (i.e., including both UV- 
bright and faint galaxies). Integrating the stellar mass 
function given by Reddy et al. (2009) suggests that while 
galaxies in the stellar mass range probed by OSIRIS (i.e., 
0.1 - 7.8 x 10 10 M Q ) constitute ~ 66% of the total space 
density of galaxies at z ~ 2, 90% of the total popula- 
tion have masses smaller than the average of the OSIRIS 
sample. 
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5.2. Expanding the Sample 

In order to construct a more complete picture of the 
overall z ~ 2 galaxy population we therefore combine 
our results with those of similar studies recently under- 
taken by a variety of authors. The largest sample of IFS 
observations of galaxies at redshift z ~ 2 to date is the 
"SINS" survey using SINFONI on the VLT, the results 
of which have been presented by Forster Schreiber et 
al. (2006); Genzel et al. (2006); Bouche et al. (2007); 
and more recently Shapiro et al. (2008) and Genzel et 
al. (2008). 15 While all of the SINS galaxies (a few of 
which are in common with the present survey) exhibit 
comparably high velocity dispersions to the 13 OSIRIS 
galaxies reported here, many of them also show ordered 
velocity shear > 100 km s _1 to which Genzel and collab- 
orators (see particularly Shapiro et al. 2008) have been 
relatively successful in fitting turbulent disk models. 

In part, the different prevalence of shear between the 
SINS and OSIRIS surveys is a consequence of their sensi- 
tivity. Due to a combination of optical design differences 
and a larger (0.1" ) spaxel scale SINFONI is roughly twice 
as sensitive as OSIRIS and therefore able to probe fainter 
surface brightness features to larger radii. Combined 
with longer integration times (18000 s. vs. 8100 s in 
the case of Q2343-BX389, for example), SINFONI is ca- 
pable of detecting large-scale velocity gradients in lower 
surface-brightness galaxies undetectable with OSIRIS. In 
contrast, the image quality of the OSIRIS data is gen- 
erally superior, with LGSAO image correction for all of 
the target galaxies and a typical PSF ~ 100 mas, com- 
pared to the large fraction of the SINS sample (13 of the 
14 galaxies discussed by Forster Schreiber et al. 2006) 
which were obtained in seeing-limited mode with a PSF 
i~ 500 mas (with some notable exceptions; e.g., Genzel 
et al. 2006). One simple explanation may therefore be 
that velocity shear is less noticeable in the inner few kpc 
of z ~ 2 star-forming galaxies (where the OSIRIS survey 
is most sensitive) and only becomes significant relative 
to the line-of-sight velocity dispersion at larger radii and 
fainter surface brightnesses (where the SINS survey is 
most sensitive). Indeed, the maximum kinematic slope 
of DSF227a-C2 (the galaxy which presents the strongest 
case for resolved velocity shear in an inclined system for 
the OSIRIS sample) ^§ ~ 35 km s" 1 kpc~ x is compa- 
rable to that of the weakest rotator in the Genzel et al. 
(2008) sample (SSA22a-MD41), and if this slope extends 
to larger radii the two galaxies may also have compa- 
rable fshear when observed to a similar limiting surface 
brightness. 

In general however, such a purely instrumental bias 
between the two surveys is an unsatisfactory explana- 
tion for the observed kinematic differences. While the 
upper end of the kinematic slopes observed with OSIRIS 
overlaps with that of the SINFONI "massive disk" popu- 
lation, the mean value for the 13 galaxies presented here 
((<r) = 19 km s" 1 kpc^ 1 ) is substantially less than the 
mean ((4r-) — 47 km s _1 kpc -1 ) of the 5 galaxies dis- 
cussed by Genzel et al. (2008). Indeed, 6 of the 13 
galaxies observed with OSIRIS have 4| consistent with 

15 We refer the reader also to Forster Schreiber et al. (2009), 
which first became available during the publication process of this 
paper. 



to within observational uncertainty. In addition, both 
SINFONI and OSIRIS find similar kinematics in the few 
galaxies (e.g. Q1623-BX502) which were successfully ob- 
served by both surveys and do not indicate substantial 
shear in a low surface-brightness component. 

Most likely, the kinematics probed by the two sur- 
veys represent a continuum ranging from genuinely 
dispersion-dominated to rotationally supported systems. 
This range in kinematic properties is perhaps unsurpris- 
ing in light of the relative physical properties of the galax- 
ies included in the two samples. As discussed above, 
OSIRIS targets were drawn from the opticaly-selected 
"BX" galaxy sample; the successful observations tend to 
have stellar masses in the less-massive to typical-mass 
range. In contrast, many of the SINS galaxies with sig- 
nificant velocity shear have inferred stellar masses that 
would place them in the top quartile of the general "BX" 
sample. Forster Schreiber et al (2006) discuss a sample of 
12 UV-selected galaxies (excluding Q1307-BM1163 due 
to its lower redshift, and Q1623-BX663 due to its sub- 
sequent identification as an AGN) with a similar range 
in M„ compared to the OSIRIS sample. The four galax- 
ies in the sample with unambiguous velocity gradients 
with v sheal > 100 km s -1 (SSA22-MD41, Q2343-BX389, 
Q2343-BX610, Q2346-BX482) have an average stellar 
mass (log(M*/M Q )) = 10.99, 8 times greater than the 
average of the rest-UV selected population. 16 In their de- 
tailed discussion of disk kinemetry, Shapiro et al. (2008) 
compare to a similarly high mass subset of 11 galax- 
ies chosen for their particularly high quality data. The 
mean stellar mass of their 5 rest-UV selected galaxies is 
(log(M*/M )) = 10.9, the remaining 6 are selected ac- 
cording to their rest-frame optical BzK color (e.g. Daddi 
et al. 2004). While there is a large amount of over- 
lap between the BX and BzK samples, the particular 
BzK objects selected all have bright K band magnitudes 
(K s < 20.5) whereas the BX objects (selected without 
regard to K magnitude) extend considerably fainter, to 
K s ~ 22 (see discussion by Reddy et al. 2006b). The 
BzK sample also tend to be slightly less strongly nu- 
cleated with fewer spatial irregularities (see discussion 
by Law et al. 2007b), consistent with the interpreta- 
tion that they might more accurately be described as 
evolved stellar "disks" than the BX sample. The Bouche 
et al. (2007) sample in turn is roughly a superset of the 
Forster Schreiber et al. (2006) rest-UV selected galaxies, 
-ftT-bright rest-optically selected galaxies, and submillimc- 
ter galaxies (e.g. Chapman et al. 2005). In turn, Genzel 
et al. (2008) focus on five primary galaxies (four pre- 
viously discussed by Forster Schreiber et al. (2006) and 
one additional galaxy BzK6004) with a mean stellar mass 
(\og(M^ / Mq)) = 11.0, 8.5 times larger than the mean of 
the rest-UV selected population, and more massive than 
> 99% of the overall z ~ 2 galaxy population (based on 
the mass function of Reddy et al. 2009) . 

These data are summarized in Table [6] and Figure 
1101 which plots the maximum observed shear velocity 
of galaxies as a function of stellar mass. While there 

16 Where available, we have used stellar masses for the SINS 
galaxies calculated by Genzel et al. 2008. Otherwise, masses are 
drawn from our own stellar population modeling. Both estimates 
are consistent and adopt a Chabrier (2003) IMF. For our galaxies 
with more than one component we focus on the velocity structure 
of the most massive component. 
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Fig. 10. — Logarithmic plot of stellar mass vs. maximum line-of-sight velocity shear for our 15 galaxies (DRL 2009) and galaxies discussed 
by Forster Schreiber et al. (2006; NFS 2006), Genzel et al. (2006; RG 2006), Genzel et al. (2008; RG 2008), and Stark et al. (2008; DPS 
2008). All data shown here are uncorrected for inclination effects. The dashed line denotes (log(<r m ean)) of our sample (note that most 
galaxies fall below this line), the dotted line represents a simple least-squares fit to the data (log (f s h C ar/km s _1 ) = 0.30x log(Af*/M0) 
- 1.39). The histogram in the lower part of the plot (arbitrary scale) represents the relative number of galaxies in each logarithmic mass 
bin for the full sample of 818 optically-selected galaxies with stellar mass information in the redshift range 1.8 < z < 2.6. Uncertainties in 
^shcar are given where available, the solid line in the lower left corner of the plot indicates the typical uncertainty in M„ . 



is considerable scatter in the diagram (unsurprisingly 
given that u s hear is n ot, and in most cases cannot, be 
corrected for inclination) we note that all galaxies with 
^shoar greater than the typical velocity dispersion have 
masses 2 — 20 times greater than the average for rest-UV 
selected galaxies at these redshifts. 

If velocity shear is indicative of systematic rotation, 
this suggests that stable rotation may be more prevalent 
(or at least more easily observed) among galaxies which 
have already accumulated a sizeable stellar population. 
It may be that the lowest mass galaxies tend to display 
negligible velocity shear and extremely high Ha surface 
brightness, while star formation in more massive galaxies 
with appreciable old stellar populations tends to trace 
the kinematics of any underlying disks with greater fi- 
delity. While this trend might hold in a statistical sense, 
however, we stress that individual galaxies do not al- 



ways fit this relation: one of our highest mass objects 
(Q2343-BX453) is perhaps the single best example of a 
galaxy without any resolved large-scale kinematic struc- 
ture, while some low mass objects (Q1700-BX490) show 
obvious structure. Indeed, the observation by Stark et 
al. (2008) of a gravitationally lensed galaxy offered an 
exceptionally fine-scale (resolution ~ 100 pc) view of a 
low-mass galaxy whose stellar mass and velocity struc- 
ture are remarkably similar to that of Q1700-BX490 (Ta- 
ble Ej) with smoothly varying velocity shear. 

5.3. Local Analogues and Galaxy-Galaxy Mergers 

Recent studies of the clustering properties of z ~ 2 
star-forming galaxies in conjunction with large cosmo- 
logical N-body simulations (e.g. Springel et al. 2005) 
have indicated (Conroy et al. 2008; Genel et al. 2008) 
that these galaxies evolve into a variety of galaxy types 
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by the present day, including typical L* galaxies. Local 
galaxies with similar star formation properties, however, 
tend to be gas-rich mergers with multiple nuclei and tidal 
features (e.g., Sanders et al. 1988; Bushouse et al. 2002), 
such as the ULIRGS Mrk 273 and IRAS 15250+3609. 
Similar to the high-redshift sample, these galaxies show 
large velocity dispersions with little resolved kinematic 
substructure in the regions of brightest Ha emission, de- 
spite relatively strong (~ 200 km s _1 ) and highly dis- 
turbed features in lower surface brightness regions (Col- 
ina et al. 2005). Perhaps the best analog of these high- 
redshift galaxies that has been found to date however 
is the supercompact subsample of UV-luminous galax- 
ies (ScUVLGs) discovered using GALEX (Heckman et 
al. 2005). These systems have a similar range of specific 
SFRs, metallicities, and dust content to high-redshift Ly- 
man Break Galaxies (see discussion by Heckman et al. 
2005; Hoopes et al. 2007; Basu-Zych et al. 2007), and re- 
cent evidence (Basu-Zych et al. 2009) indicates that their 
kinematics are also similar to those of z ~ 2 — 3 galax- 
ies With high (T m ean > 80 km S _1 and Ushcar/Cmcan < 1. 

Many of these ScUVLGs show tidal features and other in- 
dications of active mergers on small physical scales and at 
low surface brightnesses which would not be detectable at 
high redshift, leading Overzier et al. (2008) to conclude 
that rapidly star forming galaxies such as our rest-UV 
selected sample are mergers of gas-rich galaxies which 
have triggered super starbursts on scales ~ 100 - 300 pc. 

Numerous morphological studies (e.g., Conselice et al. 
2003; and references therein) have interpreted the ir- 
regular morphologies of z ~ 2 galaxies as evidence for 
galaxy-galaxy mergers in the high-redshift universe, and 
our own kinematic data may further bolster this interpre- 
tation in some cases (e.g., HDF-BX1564, Q1700-BX490). 
One natural suggestion therefore is that the observed 
multiple-clump morphologies and unusual kinematics of 
these galaxies may be a consequence of merger activity 
which prevents the formation of a stable disk at early 
times (e.g., Ostriker 1990). Such lumps can disrupt the 
preferred kinematic axis of disk galaxies in the local uni- 
verse (e.g., Benson et al. 2004), and their effects would 
likely be stronger at higher redshifts when dark matter 
haloes were merging more rapidly (e.g., Zentner & Bul- 
lock 2003; Fakhouri & Ma 2008). By inflating the ve- 
locity dispersion of putative disks in all directions (and 
not simply perpendicular to the disk) , such mergers with 
either luminous or dark satellites may help mask weak ro- 
tational signatures while more massive galaxies with cor- 
respondingly larger rotational velocities may naturally be 
more resilient to such disruption. 

However, comparison with cosmological merger trees 
suggests (Conroy et al. 2008; Genel et al. 2008) that 
major mergers are too rare to account for the high space 
density of all z ~ 2 star-forming galaxies. In addition, 
some of our best merger candidates are simultaneously 
those with the clearest rotation curves (Q1700-BX490). 
It is therefore unlikely that the uniformly high velocity 
dispersion of every single star-forming galaxy observed 
to date at redshift z ~ 2 is due to large-scale merger 
activity, suggesting that the star-formation may instead 
be due to a different triggering mechanism than observed 
in the local universe. 

5.4. An Evolving Picture of Galaxy Formation 



According to classic theories of galaxy formation (e.g., 
White & Rees 1978; Mo, Mao, & White 1998; see also 
the comprehensive review by Baugh 2006), hot-mode 
spherical accretion dominates the gas accretion history 
of galaxies. Once a sufficiently massive dark matter halo 
has virialized, gas collapses through the virial radius of 
the potential well, heating via shocks to the virial tem- 
perature of the host halo. As this gas cools over time 
(largely by collisionally excited line radiation for haloes 
with 10 4 K < T < 10 6 K, and bremsstrahlung radia- 
tion for haloes with T ~ 10 7 K; White 1994) it collapses 
to form a rotating disk supported by angular momen- 
tum which the cooled gas has been unable to shed. This 
gaseous disk gradually grows over time as gas at pro- 
gressively greater radii is able to cool and collapse, and 
is posited to be the home of the bulk of active star for- 
mation. 

The rotationally supported gas disks predicted by such 
a classical model are clearly at odds with the randomly 
directed kinematics observed in a substantial fraction of 
galaxies at z ~ 2, prompting us to revisit our under- 
standing of the complex role of gas accretion in galaxy 
formation. Unlike most galaxies in the local universe, the 
star-forming galaxies at z ~ 2 contain extremely large 
quantities of gas as compared to their stellar mass, with 
the most extreme examples inferred to have gas masses 
> 10 times larger than M* (Erb et al 2006c), and com- 
parable to their kinematically-derived dynamical mass. 
Among the UV-selected galaxies in the Ha sample of 
Erb et al. (2006c), the inferred gas fraction fi ( §3.2| de- 
creases with increasing M* ; Figure QT] recasts the trend 
of shear velocity in terms of fi. Galaxies with the high- 
est gas fraction fi > 0.5 tend to have negligible evidence 
for rotating kinematic structure and are dominated by 
randomly oriented velocity dispersions, while those with 
the lowest gas fractions fj, < 0.5 have a more obvious 
preferred axis for their angular momentum with shear 
velocities v s hear <; 100 km s _1 . We note, however, that 
while we reproduce the stellar mass - metallicity rela- 
tion of Erb et al. (2006a; Fig. Q21 left-hand panel) there 
is no clear relationship between oxygen abundance and 
apparent rotational velocity (Fig. [T2J right-hand panel) . 

Given that these galaxies have obviously amassed large 
quantities of cold gas within a very small volume (typi- 
cally r < 1 — 2 kpc) , the angular momentum of this gas 
must necessarily be small. There must exist, therefore, 
some mechanism that can funnel gas to small regions 
within typical z ~ 2 star-forming galaxies efficiently. 

Assuming for a moment that a gaseous disk is the 
first baryonic structure to form in these haloes, grav- 
itational instabilities (e.g. Ostriker & Peebles 1973) 
within this gas-dominated disk may suffice to produce 
such concentrated gas reservoirs. Indeed, Noguchi et al. 
(1999) showed that efficient cooling mechanisms could 
lead to fragmentation of gas-dominated disks into self- 
gravitating clouds with mass ~ 1O 9 M , the distribution 
and relative surface brightness of which might naturally 
give rise to an irregular, multiple component morphol- 
ogy similar to those observed for typical star forming 
galaxies at z ~ 2 (e.g., Abraham et al. 1996; Conselice 
et al. 2005; Elmegreen et al. 2005; Lotz et al. 2006; 
Law et al. 2007b). Likewise, simulations by Bournaud 
et al. (2007) also suggest that "clump-cluster" galax- 
ies (i.e. those with multiple-component morphologies) 
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Fig. 11. — Plot of maximum line-of-sight velocity shear versus gas fraction fi. Symbols are the same as in Figure HOI the dashed line 
denotes (log(<T mca n)} of our sample. 



are consistent with resulting from the fragmentation of 
unstable primordial disks. These simulated clumps typ- 
ically have gas fractions (~ 50%), star formation rates 
(~ 30 Mq yr -1 ), and lifetimes (~ 400 Myr) comparable 
to those of the z ~ 2 — 3 star forming galaxy population 
(Erb et al. 2006c). 

In gas-rich galaxies such as those of our target sample, 
these massive gaseous clumps may actually drive the dy- 
namical evolution of the entire baryonic component. As 
discussed by Bournaud et al. (2007), simulated galaxies 
during the "clump-cluster" phase are highly disturbed 
with large local velocity dispersions ~ 50 — 80 km s _1 (i.e. 
similar to those observed in our target galaxies) and neg- 
ligible coherent rotational signatures on scales ~ 2 — 3 
kpc. At later times (~ 0.5 - 1 Gyr) these clumps are 
expected to gradually disperse and sink to the center of 
the galaxy through dynamical friction, providing a possi- 
ble mechanism for bulge formation aided by the massive 
starburst (SFR ~ 100 Mq yr" 1 ) predicted to occur when 
the remaining clouds merge in the galactic center (e.g., 



Immeli et al. 2004ab; Elmegreen 2008). At such late 
times it is likely that a significant stellar population has 
been formed which can help stabilize the system against 
further instability and permit the observation of more 
regular kinematic structures (as suggested by Fig. [T0|) . 

Another possibility however is that the classical gas- 
disk phase may be bypassed entirely by cold flows of low 
angular-momentum gas accreted directly from cosmolog- 
ical filaments (e.g., Keres et al. 2005; Dekel & Birnboim 
et al. 2006; Birnboim et al. 2007). Such flows could 
be responsible for rapidly delivering large quantities of 
cold gas to the galaxy's central regions, producing phys- 
ical conditions that in the local universe are achieved 
only in gas-rich mergers. The detailed physical struc- 
ture which might result from such an accretion model is 
unclear however, especially given the limited ability of 
numerical models to resolve the relevant spatial scales. 

Regardless of the mechanism by which this low angu- 
lar momentum cold gas is acquired, our observational 
data suggest the possibility that many young z ~ 2 star- 
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the dashed line indicates solar metallicity (Asplund et al. 2004). Upper limits to the metallicity of galaxies undetected in [Nirjare indicated 
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forming galaxies could form roughly spherical early stel- 
lar populations in a rapidly-accreted pool of low angular 
momentum gas experiencing a rapidly changing gravi- 
tational potential dominated by cooling gas. As the 
gas fraction of these galaxies decreases this early stel- 
lar population may help stabilize an extended gaseous 
disk which forms at later times from the gradual, more 
adiabatic accretion of gas from the galactic halo with in- 
creasingly greater angular momentum. The dynamics of 
more massive, evolved galaxies might therefore tend, as 
we observe, to show greater rotational structure with the 
formation of such an extended gaseous disk. 

Such an interpretation may be bolstered by the obser- 
vation of a correlation between D s hear and the kinematics 
of interstellar absorption lines traced by rest-UV spec- 
troscopy (Steidel et al. in prep.). In Figure [T3l we plot 
the offset Ai^jsM-ncb) between the absorption and sys- 
temic redshifts as a function of f shear for our galaxy sam- 
ple and that of Forster Schreiber et al. (2006) for which 
we have rest-UV spectra. While the majority of z ~ 2 
galaxies have rest-UV absorption features with centroids 
blueshifted relative to the systemic velocity by ~ 200 
km s" 1 (e.g. Pettini et al. 2002; Shapley et al. 2003; 
Steidel et al. in prep.), there is a tendency for galaxies 
with the greatest observed velocity shear to have lower 
AwiSM— neb)i indicating increased optical depth to UV 
photons at low systemic velocities (as noted previously 
by Erb et al. 2006a and Law et al. 2007a). Although 
the properties of the outflowing gas traced by interstellar 
absorption lines are generally unrelated to the kinemat- 
ics of the nebular emission lines (as any nebular emission 
within the outflowing gas is too low surface brightness 
to be detectable) it is possible that galaxies with a more 
well-defined kinematic axis could contain a greater frac- 
tion of their ISM at close to the systemic redshift (see 
discussion by Steidel et al. in prep). 

This relationship between gas content and i> s hear may 
begin to break down however at lower redshift. Obser- 
vations by Wright et al. (2009) at z ~ 1.6 find that 
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Fig. 13. — Centroid of ISM absorption lines relative to the 
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some galaxies in a similar range of stellar masses and gas 
fractions to our own sample (and selected in a similar 
manner) exhibit significant kinematic gradients greater 
than their apparent velocity dispersion. The applicabil- 
ity of this relation at different redshifts and halo masses 
is therefore uncertain, but might naturally be expected 
to evolve along with changes in the cooling time of gas 
in galactic haloes. 

6. SUMMARY 

We have presented Keck-LGSAO observations of 13 
star-forming galaxies at redshift z ~ 2 — 3 which are rep- 
resentative of the typical stellar mass of the optically- 
selected galaxy sample that dominates the global star 
formation activity at these epochs. This population 
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km s that cannot be explained by beam-smearing ef- 
fects and is an intrinsic property of the gas. These ran- 
dom motions are larger than the amplitude of any sys- 
temic velocity shear w s hcar within the central 2-3 kpc, al- 
though in some cases v s hcar may increase if the observed 
velocity gradient persists to larger radii and fainter sur- 
face brightness. 

The dominance of random motions (or the absence 
of significant shear velocity) appears to be related to 
the stellar mass and inferred gas fraction: galaxies 
with larger velocity gradients tend to have larger stel- 
lar masses (2 — 20 times larger than the mean for the 
parent sample), with smaller values of the gas fraction 
(p < 0.5). While there is no similar relation between 
metallicity and kinematics, there is some evidence to sug- 
gest that galaxies with the greatest large-scale velocity 
structure may contain a greater fraction of their total 
interstellar medium close to the systemic redshift. 

Although the observed kinematics of the high redshift 
galaxies are similar to those of merger-driven starburst 
galaxies in the local universe, it is unlikely that merg- 
ers are the primary reason for the observed dispersion- 
dominated kinematics at z <~ 2 — 3. Instead, more effi- 
cient cooling and settling of larger amounts of low angu- 
lar momentum gas (whether attributed to "cold accre- 
tion" or some other process) may generically give high 
redshift galaxies properties that arc produced only in 
mergers in the local universe. As the large, centrally con- 
centrated gas supply is rapidly converted into stars, the 
growing stellar population may help stabilize the galaxy 
and permit the formation of stable disk-like structures 
from higher angular momentum gas accreted over longer 
timescales. 

While such a qualitative picture remains speculative 
given the current state of both observations and theory, 
it is without doubt an over-simplification to classify high 
redshift galaxies as either "disks" or "mergers" ; in fact, 
many, and possibly even most, arc neither. High gas- 
phase velocity dispersions (and small velocity gradients) 
are evidently a natural consequence of the instablitics 
resulting when cold gas becomes dynamically dominant, 
as it seems to be in the central few kpc of a large fraction 
of star- forming galaxies at z ~ 2 — 3. Learning more 
precisely how and why this occurs may be the key to 
understanding galaxy formation. 
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QUANTIFYING THE MEAN PROPERTIES OF 
INCLINED SYSTEMS 



In cases for which there is obvious morphological or 
kinematic symmetry indicative of an inclined system, 
foreshortening of the respective contours may be used 
to correct apparent rotational velocities for the inclina- 
tion with confidence. In other cases however the degree 
of such foreshortening may not be obvious, and it is fre- 
quently popular to apply an 'average' inclination correc- 
tion to determine the intrinsic rotational velocity in a 
statistical sense. Confusingly, the value of this average 
correction frequently differs in the literature, with some 
authors quoting (i) — 45° and (sinz) = f = 0.64 (e.g., 
Erb et al. 2006c; Bouche et al. 2007) and others adopt- 
ing (i) = 57.3° and (sini) = f = 0.79 (e.g. Rix et al. 
1997; Weiner et al. 2006). In an effort to foster consis- 
tency we present a simple derivation of the correct factor 
instead of simply stating an adopted value. 

The easiest way of visualizing the problem is to imagine 
the range of possible inclinations i at which a disk can be 
viewed. While i ranges from 0° to 90° , the average value 
of i is not 45°, since there are many more ways in which 
a disk may be viewed from its edge than from its pole. 
Explicitly, we frame the problem as follows. Given a 
collection of disks oriented isotropically in space, we wish 
to determine the mean expected inclination i between 
our line of sight to a given disk and the vector normal 
to the disk (where i — 0° represents a disk viewed face- 
on) . It is simplest to reduce this problem to the case of a 
single disk and consider it from the reference frame of the 
disk itself. Assuming that a viewing location is chosen 
at random, the distribution of possible lines of sight will 
uniformly cover the sky as seen from the reference frame 
of the disk. The probability dP that a viewer is located 
within a given patch of sky is therefore proportional to 
the differential solid angle d£l subtended by that patch 
of sky. Adopting a spherical polar coordinate system 8, (j> 
centered on the disk where 4> is the azimuthal angle and 8 
(ranging from — ^ to \ ) is the polar angle the differential 
solid angle is given by 



dfl = cos 8 d8 d(j). 

It is therefore possible to calculate the 
value (i) by integrating over the entire sky 
symmetry present in the problem, we may 
tegral over the azimuthal coordinate <p and 
polar integral to consider only the range 8 = 
ducing the differential solid angle to dfl = 
this range, the inclination is related to the 
by 8 = j — i. Therefore, 



(3) 

expectation 
Noting the 
omit the in- 
collapse the 
= to § , re- 
cos 8d8. In 
polar angle 



J 2 i cos 8 d8 
/ f cos 8 d8 ' 



(4) 



Substituting the identities cos 6* = sini and d6 = —di 
we obtain 

= tfi sin i di = [sini -i cos i] | = ^ y 
j 2 sin i di [— cos i] 2 

Since the observed velocities of an inclined disk are ef- 
fectively "foreshortened" by a factor sin i, it is usually of 
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interest to calculate the mean velocity reduction factor. 
Similarly to the calculations performed above, we find 

L 2 sin i cos 9 dO f„ 2 sin 2 idi n 
(sinj) = Jo_ = Jo_ = 7 = 0.79. 

J Q 2 cos 9 d9 J 2 sin i di 4 

(6) 
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2700 


Kn2 


Ha 


50 


70 




















2008 


Jun 


7200 e 


Kn2 


Ha 


100 


80 


01623-6X455' 


2.4079 


16 


25 


51 


664 


+26 


46:54 


60 


2008 


Jun 


5400 


Kn4 


Ha 


50 


85 




















2007 


Sep 


3600 


Hn4 


[Om] 


50 


70 


Q1623-BX663 


2.4333 f 


16 


26 


04 


586 


+26 


47:59 


80 


2007 


Jun 


5400 c 


Kn4 


Ha 


50 


100 


Q1700-BX563 


2.292S 


17 


01 


15 


875 


+64 


10:26 


15 


2007 


Jun 


1800 c 


Kn3 


Ha 


50 


210 


Q1700-BX691 


2.1895 f 


17 


01 


06 


117 


+64 


12:09 


70 


2006 


Jun 


10800 h 


Kn2 


Ha 


50 






















2007 


Jun 


5400 c 


Kn2 


Ha 


50 


200 


Q2206-BX102 


2.2104 f 


22 


08 


50 


751 


-19 


44:08 


24 


2007 


Sep 


5400 


Kn2 


Ha 


50 


90 


Q2343-BX389 


2.1716*' 


23 


46 


28 


911 


+12 


47:33 


90 


2007 


Jun 


8100 e 


Kn2 


Ha 


50 


160 


Q2343-BX442 


2.1760* 


23 


46 


19 


362 


+ 12 


48:00 


10 


2007 


Sep 


9000 c 


Kn2 


Ha 


50 


140 


Q2343-BX587 1 


2.2429 


23 


46 


29 


192 


+12 


49:03 


71 


2006 


Oct 


5400 


Kn3 


Ha 


50 


105 



a Vacuum heliocentric rcdshift of primary nebular emission line. b Total observing time, mean value for detected sources was ~ 2 
hrs. c Primary targeted emission lino. d FWHM of the AT-band PSF (mas) during on-axis TT star observation (before/after spatial 
smoothing respectively). Poor observing conditions/ Rcdshifts estimated from NIRSPEC spectra. g Redshift estimated from rest-UV 
spectrum. 11 Individual exposures were each 300 seconds. 1 Galaxy detected, but quality too poor for analysis. 



28 



KINEMATICS OF HIGH REDSHIFT GALAXIES 



TABLE 2 
Nebular Line Fluxes 





Aneb a 








[O III] 


L e 




[N II] 


[N II) 




Galaxy 


(A) 


2 neb 




Fdke" 


(A5007) 


(10 42 erg s 1 ) 


(A4960) 


(A6549) 


(A6585) 


12+ \og(0/Hy 


OfUzLQ RYcn 
v^U4:4:y-DAyo 


19737.8 


2.0067 


6.8 ± 0.2 






3.0 ± 0.1 




1.3 ± 0.4 


0.8 ± 0.3 


o.o t in u.uy 


Q1217-BX95 


22479.6 


2.4244 


6.5 ± 0.4 






6.4 ±0.4 




< 1.2 


< 1.2 


< 8.48 


HDF-BX1564 


21156.7 


2.2228 


9.7 ±0.8 


17.2 ± 1.4 




7.0 ±0.4 




< 1.5 


3.9 ±0.5 


8.67 ± 0.04 


Q1623-BX453 


20888.3 


2.1820 


16.4 ±0.4 


27.6 ±0.4 




12.8 ±0.3 




3.0 ±0.4 


6.5 ±0.4 


8.67 ± 0.02 


Q1623-BX502 


20715.8 


2.1557 


10.0 ±0.2 


26.4 ±0.8 




5.3 ±0.1 




< 0.7 


< 0.7 


< 8.24 


Q1623-BX543 


17634.4 


2.5211 




17.2 ± 1.4 


25.0 ±0.5 


42.5 ±0.8 


8.0 ±0.5 








Q1700-BX490 


22291.8 


2.3958 


31.6 ±0.8 


35.4 ± 1.2 




34.6 ± 0.9 




< 2.4 


< 2.4 


< 8.26 


Q1700-BX710 


21628.4 


2.2947 


4.5 ±0.2 






3.4 ±0.2 




< 0.6 


0.6 ±0.2 


8.40 ±0.08 


Q1700-BX763 


21609.8 


2.2919 


2.2 ±0.3 






1.3 ±0.2 




< 0.9 


0.3 ±0.3 


8.41 ±0.25 


DSF2237a-C2 


21621.8 


3.3172 






7.9 ±0.4 


16.8 ±0.9 


2.0 ±0.3 








Q2343-BX418 


21698.2 


2.3053 


8.2 ±0.4 


16.0 ±0.4 




3.8 ±0.2 




< 1.2 


< 1.2 


< 8.42 


Q2343-BX513 


20404.3 


2.1082 


5.4 ± 0.2 


20.2 ±0.8 




2.9 ±0.1 




< 0.6 


1.8 ±0.2 


8.63 ±0.03 


Q2343-BX660 


20835.4 


2.1739 


11.2 ±0.2 


18.8 ±0.8 




4.2 ±0.1 




< 0.6 


< 0.6 


< 8.18 



a Vacuum heliocentric wavelength of primary nebular emission: [Om] A5007 for DSF2237a-C2 and Q2343-BX415, [Nil] for Q1623-BX455, Ha for all 
others.^ 3 Heliocentric rcdshift of primary nebular emission linc. c Emission line flux in units of 10 — 17 erg s _1 cm -2 . Uncertainties quoted arc la and 
based on random errors, global systematic uncertainty is ~ 30%. Limits represent 3 a limits. 1 ^ Ha emission line flux observed by Erb et al. 2006c 
(includes a factor of 2 correction for aperture losses). Extinction-corrected primary nebular emission line luminosity (Ha or [Olll]). Uncertainties 
quoted are based on random errors, global systematic uncertainty is ~ 30%. f Oxygen abundance using the N2 calibration from Pettini & Pagel (2004). 
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TABLE 3 
OSIRIS Morphologies 







r b 


d 2c c 








Galaxy 


(kpc^) 


(kpc) 


(kpc) 


G 




M20 


O0449-BX93 


5.0 ± 0.8 


1.3 ± 0.1 




0.24 


1.2 


— 1.08 


Ol 21 7-RXQ 1 ! 


1.0 ± 0.4 


0.6 ±0.1 




0.13 


0.8 


—0.89 


HDF RX1 564 


< 3 (SE) 


1 q j- n 9 CNW1 
< 1 (SE) 


7.0 


0.10 


13.6 


—0.51 


Q1623-BX453 


6.1 ± 0.5 


Ta ± 0.1 




0.22 


1.3 


-1.36 


Q1623-BX502 


6.2 ± 1.3 


1.4 ±0.2 




0.21 


0.7 


-1.48 


Q1623-BX543 


3.6 ±0.5 (N) 
1.5 ±0.5 (S) 


1.1 ±0.1 (N) 
0.7 ±0.1 (S) 


6.7 


0.19 


9.7 


-1.09 


Q1700-BX490 


8.3 ±0.4 (W) 
1.6 ±0.4 (E) 


1.6 ±0.1 (W) 
0.7 ±0.1 (E) 


3.4 


0.15 


5.8 


-1.40 


Q1700-BX710 


4.1 ±0.5 


1.1 ±0.1 




0.23 


0.5 


-1.51 


Q1700-BX763 


3.9 ±0.9 


1.1 ±0.1 




0.14 


1.5 


-1.22 


DSF2237a-C2 


2.7 ± 0.4 


0.9 ± 0.1 




0.17 


2.4 


-1.16 


Q2343-BX418 


2.2 ±0.6 


0.8 ±0.1 




0.13 


1.8 


-1.32 


Q2343-BX513 


6.2 ±0.7 


1.4 ±0.1 




0.19 


1.1 


-1.28 


Q2343-BX660 


8.5 ±0.5 


1.6 ±0.1 




0.18 


1.7 


-0.99 



a Area of nebular emission. Uncertainty represents half the PSF correction. Individual 
components are identified in brackets.* 5 Radius of nebular emission. Uncertainty repre- 
sents half the PSF correction. Individual components are identified in brackets. c Dis- 
tance between the two morphological components. d Gini. e Multiplicity. 
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TABLE 4 
Kinematic Properties 



Galaxy 


rr a 

V mean 

(km s _1 ) 


Onet 

(km s _1 ) 


n , c 

'-'shear 
(km S _1 ) 


V'2c 

(km s _1 ) 


A c 
^(ISM-neb) 

(km s _1 ) 


I : \ 

( ) 


^shear / "mean 


-™dyn 
tt r»10 Ti/T \ 

(1O 1U M ) 


OfUzLQ RYQ3 


56 ± 20 


71 ± 5 


13 ± 8 




—270 


330 


0.2 i 0.2 


0.5 




69 ± 23 


61 ± 5 


12 ± 13 




—298 


110 


0.2 ± 0.2 


0.3 


HDF-RX1 564 


fin + 1 

UU _1_ -Lijll'l VV J 


103 ± 14 


19 + 12 fNWl 


171 


—447 


325 CNW1 

*j£i\J 1 11 VV / 


2 + 02 (NWl 


5 CNW1 

yj •<j y 1 » v v / 




ftp. 4. sfSFA 




tj _1_ luLl 1 






325 CSF/) 


U . J- _1_ \J • J. I kji— 1 J 


< 9 fNW1 


Q1623-BX453 


78 ± 23 


94 ± 5 


29 ± 10 




-849 


45 


0.4 ± 0.2 


1.0 


Q1623-BX502 


67 ±20 


68 ±3 


35 ± 11 




-67 


30 


0.5 ± 0.3 


0.7 


Q1623-BX543 


139 ± 32 (N) 


153 ± 7 


39 ± 4 (N) 


125 


-94 


315 (N) 


0.3 ±0.1 (N) 


2.5 (N) 




60 ± 10 (S) 




22 ± 9 (S) 






225 (S) 


0.4 ±0.2 (S) 


0.3 (S) 


Q1700-BX490 


89 ± 25 (W) 


122 ±6 


50 ± 10 (W) 


109 


106 


200 (W) 


0.6 ±0.2 (W) 


1.5 (W) 




65 ± 29 (E) 




22 ± 9 (E) 






270 (E) 


0.5 ± 0.3 (E) 


0.3 (E) 


Q1700-BX710 


68 ±25 


67 ±7 


15 ± 16 




118 





0.2 ± 0.3 


0.6 


Q1700-BX763 


50 ± 14 


53 ± 5 


14 ± 16 




-264 


320 


0.3 ± 0.3 


0.3 


DSF2237a-C2 


89 ±20 


99 ± 10 


51 ± 14 




125 


220 


0.6 ± 0.2 


0.8 


Q2343-BX418 


61 ± 17 


70 ±5 


23 ± 12 




-209 


70 


0.4 ±0.2 


0.3 


Q2343-BX513 


100 ± 25 


102 ±6 


65 ± 16 




-19 


180 


0.7 ±0.2 


1.6 


Q2343-BX660 


71 ± 16 


73± 


40 ±9 




-274 


45 


0.6 ± 0.2 


0.9 



a Flux-weighted mean and standard deviation of velocity dispersions in individual spaxels in the OSIRIS map. Velocity dispersion of the 
spatially collapsed object spectrum. Uncertainty quoted is based on Monte Carlo analysis of synthetic spectra. c Shear velocity f s hear — ■^■(' t, max — 
u m i n ).^ Velocity difference between two discrete components. 6 Kinematic offset between ISM absorption line rcdshift and nebular emission line 
redshift. f Position angle of steepest velocity gradient. 8 Single-component dynamical mass within the radius probed by nebular emission. 
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TABLE 5 
Stellar Population Parameters 











Agc c 


SFR SED d 


SFR ne b° 


SFR DKE f ' 




M gas h 




Galaxy 


K 


a 

3 


(10 10 M o ) 


(Myr) 


(M yr- 1 ) 


(M Q yr- 1 ) 


(M yr- 1 ) 


e(b - vy 


(1O 1O M ) 




Detections 


Q0449-BX93 


20.51 


20.59 


0.9 


203 


47 


13 




0.135 






Q1217-BX95 


20.69 


20.79 


3.0 


571 


52 


28 




0.235 






HDF-BX1564 


19.62 


19.70 


2.9 


571 


53 


31 


34 


0.205 


3.6 


0.55 


Q1623-BX453 


19.76 


20.00 


3.1 


404 


77 


57 


100 


0.245 


3.3 


0.52 


Q1623-BX502 


22.04 


22.69 


0.3 


203 


15 


24 


79 


0.130 


3.5 


0.92 


Q1623-BX543 


20.54 


20.66 


0.4 


10 


431 


189 


98 


0.285 


4.6 


0.92 


Q1700-BX490 


19.99 


20.26 


0.4 


10 


441 


154 


160 


0.280 


5.6 


0.93 


Q1700-BX710 


20.18 


20.23 


4.4 


1278 


34 


15 


59 


0.195 


2.5 


0.36 


Q1700-BX763 


20.94 


20.98 


1.4 


905 


15 


6 


23 


0.125 


0.8 


0.36 


Dor zzo i a-Uz 


20.53 


20.71 


2.9 


255 


112 


39 




0.175 






Q2343-BX418 


21.88 


22.29 


0.1 


102 


11 


17 


32 


0.030 


1.4 


0.93 


Q2343-BX513 


20.10 


20.16 


7.8 


3000 


26 


13 


44 


0.155 


1.8 


0.19 


Q2343-BX660 


20.98 


21.21 


0.8 


1609 


5 


19 


30 


0.005 


1.8 


0.69 










N 


on-Dctcctions 












Q0100-BX210 






0.8 


102 


74 






0.215 






HDF-BX1311 


20 


48 


0.9 


255 


35 




48 


0.110 


2.3 


0.72 


HDF-BX1439 


19 


72 


4.6 


2200 


21 




48 


0.175 


3.7 


0.45 


Q1623-BX455 


21 


56 


0.8 


45 


175 




169 


0.370 


5.4 


0.87 


Q1623-BX663 j 


19 


92 


13.2 


2000 


21 




50 


0.135 


2.7 


0.17 


Q1700-BX563 






3.3 


641 


51 






0.220 






Q1700-BX691 


20 


68 


6.0 


2750 


22 




36 


0.285 


2.1 


0.26 


Q2206-BX102 






11.8 


2750 


43 






0.335 






Q2343-BX389 


20 


18 


7.2 


2750 


26 




80 


0.265 


4.7 


0.39 


Q2343-BX442 


19 


85 


11.8 


2750 


43 




44 


0.285 


3.2 


0.21 


Q2343-BX587 


20 


12 


3.3 


719 


46 




32 


0.175 


2.1 


0.39 



a Observed Vega magnitude, values in brackets are corrected for line emission (where observed). 13 Typical uncertainty < <tm* /M* >— 0.4. c Typical 
uncertainty < crAge/Age >— 0.5.^ SFR derived from stellar population model. Typical uncertainty < (Tsppt/SFR >— 0.6. e Extinction-corrected 
nebular SFR. Values for Q1623-BX543 extrapolated from [O III]/ Nebular SFR estimate from long-slit observations by Erb et al. (2006b), corrected 
for extinction and slit losses. g Typical uncertainty < (?e(b — V)/E(B — V) >— 0.7. 11 Gas masses taken from Erb et al. 2006c. 1 Gas fraction^ AGN 
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TABLE 6 
Other Galaxy Samples 









^shcar 


Af(ISM~nob) d 








Galaxy 


z a 




(km s- 1 ) 


(km s- 1 ) 


(1O 1O M ) 




12+ log(0/i?)« 


Forstcr Schrciber et al. (2006) 


Q1623-BX376 


2.4087 


20.84 


30 


-62 


0.5 


0.81 


8.43 ±0.17 


Q1623-BX455 


2.4071 


21.56 


55 


-9 


0.8 


0.87 


8.41 ±0.12 


Q1623-BX502 


2.1555 


22.04 


50 


-48 


0.3 


0.92 


< 8.16 


Q1623-BX528 


2.2684 


19.75 


65 


-220 


12.2 


0.18 


8.62 ±0.04 


Q1623-BX599 


2.3313 


19.93 


50 


-117 


6.7 


0.37 


8.55 ±0.08 


SSA22a-MD41 


2.1710 


20.5 


140 


189 


2.7 h 


0.43 h 


8.43 ± 0.08 


Q2343-BX389 


2.1728 


20.18 


235 


-76 


7.0 h 


0.22 h 


8.54 ±0.08 


Q2343-BX610 


2.2102 


19.21 


165 


-112 


17 h 


0.15 h 


8.68 ±0.01 


Q2346-BX404 


2.0298 


20.05 


20 


-277 


40 


0.31 


8.50 ±0.07 


Q2346-BX405 


2.0308 


20.27 


30 


-79 


1.6 


0.56 


< 8.16 


Q2346-BX416 


2.2406 


20.30 


70 


37 


2.5 


0.52 


8.46 ±0.07 


Q2346-BX482 


2.2569 




100 


-175 


6.6 h 


0.31 h 


8.49 ±0.07 








Genzel et 


al. (2006) 








BzK15504 


2.3834 


19.2 


171 




7.7 


0.36' 




Genzel et al. (2008) 


SSA22a-MD41 


2.172 


20.5 


165 


189 


2.7 


0.43 h 


8.27 ±0.06 


Q2343-BX389 


2.174 


20.18 


261 


-76 


7 


0.22 h 


8.51 ±0.04 


Q2343-BX610 


2.211 


19.21 


158 


-112 


17 


0.15 h 


8.66 ±0.02 


Q2346-BX482 


2.258 




213 


-175 


6.6 


0.31 h 


8.35 ±0.06 


BzK6004 


2.387 


18.9 


146 




58 


0.05 h 


8.69 ±0.02 


Stark ct al. (2008) 


J2135-0102 


3.07 




55 




0.6 




8.6 



a Systemic rcdshift, values drawn from the indicated references. Vega magnitude, uncorrected for line emission. 
Values based on our own observations except for BzK15504, BzK6004, SSA22a-MD41. c Values based on indicated 
references: v r from Forstcr Schrciber et al. (2006), v c sini from Genzel et al. (2006) and Stark ct al. (2008), Vd sin-i 
from Genzel et al. (2008). ^ Kinematic offset between ISM absorption line rcdshift and nebular emission line rcdshift, 
calculated from our rcst-UV spectra. Stellar masses for Forstcr Schrciber ct al. (2006) galaxies from our broadband 
photometry (except where indicated), for all other samples drawn from the indicated reference* Combines tabulated 
M* with M gas from D.K. Erb (priv. comm.) except where indicated. g Oxygen abundance drawn from the indicated 
reference.* 1 Values from Genzel ct al. (2008).* Values from Genzel et al. (2006). 
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